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CHAPTER  I 
INTRODUCTION 

In  |_2lJ  Payne  and  Weinbeger  presented  a  method  for 
determining  pointwise  bounds  for  the  solution  of  Dirichlet 
and  mixed  type  boundary  value  problems  for  arbitrary  second 
order  elliptic  equations  and  rather  general  domains.   In 
the  present  paper  the  writer  will  obtain  similar  results  in 
the  parabolic  case  for  Dirichlet  and  mixed  type  boundary 
value  problems  where  the  given  data  is  a  function  of  the 
domain  for  some  problems  while  for  others  the  data  may  depend 
on  both  the  solution  and  the  domain  of  the  problem. 

The  methods  of  this  paper  yield  error  estimates  in 
terms  of  the  integral  of  the  square  of  known  functions  and 
thus  the  results  lend  themselves  to  improvement  by  means 
of  the  Rayleigh-Ritz  technique  when  the  given  data  is  a 
function  of  domain  (the  linear  case).   In  such  cases  it  is 
an  unfortunate  fact  that,  even  if  the  set  of  functions  to 
which  one  applies  the   Rayleigh-Ritz  technique  is  complete 
in  the  function  space  of  the  solution,  the  actual  solution 
of  the  problem  can  not  be  obtained  using  our  present-day 
computers  due  to  the  machine  creation  of  large  round  off 
errors.   It  is  quite  possible  however  that  future 
generations  of  these  computing  devices  will  eliminate 


such  errors  and  thus  make  the  techniques  or  this  paper  of 
greater  practical  importance.   In  this  regard  we  note  that 
the  current  electronic  computing  machines  have  made  of 
practical  interest  problem  solving  techniques  which  were 
previously  only  of  theoretical  interest  since  feasible  uses 
were  not  apparent  at  the  time  of  their  evolution. 

For  many  of  the  problems  considered  in  this  thesis, 
the  existence  of  a  solution  has  not  been  established  under 
the  general  conditions  necessary  for  the  existence  of  the 
pointwise  bound.   When  certain  addition  requirements  or 
modifications  are  placed  on  the  boundary  value  problem,  by 
referring  to  the  literature  [l]  ,  |Y|  ,  [9]  to  [l5]  $     fez]    and 
[23]  we  note  that  the  existence  of  the  solution  is  guaran- 
teed. 

As  was  pointed  out  above,  the  estimates  of  this  paper 
are  in  terms  of  the  integrals  of  known  functions.   For 
another  method  of  estimating  the  solution  function  of 
boundary  value  problems,  the  reader  is  referred  to  the 
excellent  works  of  Fichera  [6 J  to  |_8J  where  the  solution 
function  is  bounded  in  terms  of  the  maximum  norm.   Using 
the  results  of  Fichera  as  a  starting  point,  Lieberstein  |l9j 
derived  a  numerical  procedure  by  which  error  bounds  were 
obtained  using  a  digital  computer. 

In  Chapter  II  we  obtain  some  preliminary  results 
which  shall  form  the  basis  of  many  of  the  conclusions  of 
the  following  chapters.   The  identity  of  Hormander  [l6]  9 


which  was  used  by  Payne  and  Weinberger  [21]  „  is  modified  in 
order  to  allow  its  use  in  a  Dirichlet  type  parabolic  problem 
and  functions  analogous  to  (i+ » X )  and  (I4. «X8)  of  [21]  are 
introduced  and  investigated..   Chapter  III  concerns  itself 
with  non-linear  normally  parabolic  problems  of  both  the 
Dirichlet  and  mixed  types.   For  the  purposes  of  this  thesis 
a  normally  parabolic  differential  operator  has  the  form 

1.1  J(w)  =^>  XfaiJ  ij"  V  z  i2 

ij  *       b  '     by 

where  the  matrix  a-^J  is  positive  definite  and  Z  is  strictly 
positive.   We  shall  see  that  the  desired  estimates  are 
obtained  by  multiplying  the  differential  operator  by  a 
parametrix  function  and  then  using  Green's  theorem  to 
generate  the  solution  at  the  point  in  question,,   The  unknown 
terms  in  the  resulting  expression  are  estimated  as  a  con- 
sequence of  the  results  of  Chapter  II  and  by  further  use  of 
the  divergence  theorem,,   In  Chapter  IV  the  problems  considered 
are  degenerate  in  character,  thus  the  matrix  a  J  becomes 
positive  semi-definite  for  some  problems  while  the  function 
Z  need  only  be  non-negative  for  others.   The  next  to  the  last 
chapter  of  the  paper  outlines  a  method  by  which  solution 
bounds  may  be  generated  for  a  class  of  elliptic  problems 
which  were  not  considered  by  Payne  and  Weinberger  in  their 
investigation  of  such  problems  [2lJ  .   In  Chapter  VI  the 
results  of  the  thesis  are  applied  to  obtain  bounds  for  the 
solution  of  a  specific  boundary  value  problem., 


CHAPTER  II 

PRELIMINARY  RESULTS 

1.   General  Definitions 

Let  V  be  an  open  (N  +  1)  dimensional  domain  of  real 
variables  (x,  y)  =  (x  ,  x2  •  .  .  xN,  y)  which  is  bounded  by 
two  hyperplanes,  y  =  o  and  y  =  7o>0,  and  a  surface  S  lying 
between  these  planes o   As  usual  we  denote  the  closure  of  V 
by  V  and  assume  that  the  divergence  theorem  holds  in  V.   The 
reader  is  referred  to  Chapter  IV  of  |_l8J  for  information 
relating  to  the  most  general  regions  in  which  this  theorem 
is  valid.   Let 

2.1  Dty)  =  {(x,y)  y  -«j ,  (x,y)ev  -  s) 
then 

2.2  V  =  VUD(o)UD(yQ)US  =  VUD(o)UD(y0)US  . 

On  the  boundary  of  V  the  Euclidean  outward  normal  is  denoted 
by 

2.3  (ni5n   )    =   (n^,    n2   .    .    .   nN,    n  ) 

where 

o      N    o 
2.k  nd  +S"nf  =  1  , 

y    V  i 


Since  V  is  a  region  in  which  the  divergence  theorem  is  valid, 
it  is  clear  that  2.3  is  (uniquely)  defined  almost  everywhere 
on  the  boundary  of  V.   In  addition,  unless  otherwise  indicated, 
we  require  that  S  be  such  that 

N 

^  2 
2.5         min^>_n,  -  m^>o 

S  1 

for  all  (x,y)6S  at  which  the  normal  vector  is  defined. 

A  function  2  -   Z(x  y)  is  said  to  be  piecewi.se  continuous 
in  V  if  V  may  be  divided  into  a  finite  number  of  subregions 
in  each  of  which  the  function  is  continuous  and  has  a  finite 
limit  as  the  boundary  is  approached.   The  above  is  the  usual 
definition  of  piecewise  continuity;  the  following  non- 
standard definitions  shall  also  be  employed  since  they 
greatly  simplify  the  presentation  of  the  problem.   When  a 
function  2(x,y)  is  said  to  be  continuously  (piecewise  con- 
tinuously) differentiable  in  x  for  (x,y)6V  then  Z  is  con- 
tinuous on  D(y)  for  o<y<y0  and  is  piecewise  continuous  in  V 

No 

while  the  $=-±9    I  =  1,  2  .  .  .  N,  are  continuous  (piecewise 

oX 

continuous)  on  D(y)  for  o<y<y  and  piecewise  continuous  in 

V.   Thus  if  Z  is  twice  continuously  differentiable  in  x  and 

piecewise  continuously  differentiable  in  y  for  (x?y)eV  then 

Z  is  continuous  in  V,  the  2JL ,  i  =  1,  2  .  .  .  N  are  con- 

^x1 

tinuous  on  D(y)  for  o<y<y  and  piecewise  continuous  in  Vs 

the  ^.   •  ,  ij  =  1,  2  o  .  .  N,  are  continuous  on  D(y) 

oxW  . 

for  o<y<y_  and  piecewise  continuous  in  Vs  and  finally  s_ 
o  ()  y 


is  piecewise  continuous  over  the  entire  volume. 

2.   The  Parametrix  of  the  Parabolic  Problem 

For  (x,y)€V  let  the  following  differential  operator  be 
defined 

2o6         im    =  (a^^j)^  +  Sy  (Z?) 

where  the  symbol  ,^  indicates  partial  differentiation  with 
respect  to  x  ,  the  summation  convention  is  used  throughout 
and  the  components  of  the  symmetric  matrix  a  J  =  a 3(x,y) 
are  piecewise  continuously  differentiable  in  x   for  (x*y)£V 
and  are  such  that 

2.7       <><«of  Si**1J5i?j«*i   fit 

for  any  numbere  (§-,,  .  .  • J  w)  and  all  (x,y)6  V  with  equality 
on  the  left  holding  if  and  only  if  }j  =  o  for  i=l,  2  .  .  .  N. 
The  function  Z  =  Z(x,y)  is  piecewise  continuously  differentiable 
in  y  for  (x,y)€V.   For  the  point  p  =  (x0,yoKD(yo)  ,  at  which 
a  bound  for  the  solution  of  certain  boundary  value  problems 
shall  later  be  obtained,  we  assume  the  functions  a1^  and  Z 
satisfy  a  Lipschitz  condition  for  (x,y)6V;  that  is,  we 
suppose  there  exist  positive  numbers  A  •  and  A  such  that 


2.8 


iij°(x,y)  -  a^(p)|  {(yQ-y)2  +  X  (x^xj)2}  ^A1^ 


2.9 


z(x,y)-  z(p)  |{(y0-y)2  ♦  lEu1-**)}  4a 


for(x,y)ev. 


The  parametrix  function  ft      is   defined  with  respect 


to   the  point  p  by 


2.10 


for 


where 


P    L 


**TT  (y    -y) 
o 


-i  U 
J 


Z(p)    2    \/£Tp5    exp[-§Up)p  (i+Cy^y)}"1 


2.11  p  =   (xoS   y0)^D(y0),    (xsy)£V,    and  Z(p)>o 


and  a(p)  is  the  determinant  of  the  inverse  matrix  a»  .(p). 
On  the  surface  D(y  )  -  Vp)  9    the  parametrix  and  its  deriva- 
tives are  defined  as  the  limits  of  these  functions  as 
(x,y)-»(x,y0)  with  (x5y)6V  and  x^0o 

Theorem  2.L   The  parametrix  6p  has  the  following 
properties  2 

a)  At  all  points  (x5y)£V  -  {p1  »  ^ne  parametrix  is 
twice  continuously  differentiable» 

b)  The  function 


2.13    (y0-y) 


4X 


>*p  [$  (y0-yf  *  2(p),o2  j(*  ) 


is  integrable  over  V  for  5  anf3  A.  satisfying 


2. 


1^ 


5<i 


and 


o  <  A  . 


8 


c)   Pop  any  bounded  function  r  which  is  continuous 
on  Vl/D(y0)  the  parametrix  is  such  that 

2-15  U.  l/J2^  -  m-  nK°^- 

Proofs  Prom  the  definition  of  ^p,  it  is  obvious  that  part 
(a)  is  satisfied.  To  establish  part  (b)  of  the  theorem  we 
note  that  for  c . 5 i  -  2S  3  .  .  .  appropriately  defined  (the 
constant  c^  is  reserved  for  later  use) 

2.16      J/i^-^  +  ^[%-^'?^],j(^)aiv 


-  If  4^ V  #<VW>1  *p[-(&  "  9  «f^  V  W)^*.  ^ 


where 


2ol? 


^(i-S)*(f)^^o-#]^«t<j^^ 


N 


r2  -  X  C*1   "  *X) 


2 


2.18  g.1J(xk-xk)+6iJ(y  -7)-a13(x,y)-a1J(p)      1,  J=l,2   .    .    .   N 

k  o       y       o 


with 


2.19 


and 


*i*  -- e*\* )[%-i?+i^-*'f\Vk%)-*%i    -a**/a a/ 


^^V^V^^'^T'^M)-^)}      *,$»/, a/ 


2.20        h,  (x^x1)  +  h  (yrt-y)  =  S(xsy)   -  «(p) 


0 


with 


21       { 


(V  <*w  Wa +&Hef]~  W  -zcf>i 


*a', v 


4=  (^[(^B^J  {e^)-*fP)} 


^ 


10 

Now  since  the  a1^  and  Z  satisfy  2.8  and  2-9>  the  functions 

gr^s  fc*^»  h.  and  h  are  bounded  over  V  and  thus  the  constants 
°k    y    i      y 

on  the  right  of  2»l6  exist.   By  integration  by  parts  with 
respect  to  r  and  then  j9    it  is  easy  to  establish  (using  2.lij.) 
that  each  of  the  integrals  on  the  right  of  2.l6  is  bounded 
and  thus  part  (b)  of  the  theorem  is  proved. 

To  establish  the  validity  of  2.l£  we  suppose  that 
£Q> o  is  chosen  so  that  for  o<£<£  we  have 


2,22 


r  rL 
=  &**      J    J    ViX^  ds  . 


Since  2.22  is  true  for  every  £  such  that  o<«£<£  and  since 
Z  and  f  are  continuous  at  p6D(y0),  we  have  that 

2.23     i^    If  f?^i^ 

£-*o   ^"^   <4/  o       I 


11 


v-a. 


UJ      0 
n 


to 


1H 


V^     Hi 


+  &*jA^o  I   I  l<,N-l)ry    C^r^Y^^^'^^Y^r^j^J.uj 


Hence   if  N  is   even  we  obtain 


2-^     iL,  n  *ix»^5 


12 


l-*o       %*fy     %, 


^u-Kf>)^\^o-^: 


cL^s  cL 


u) 


-£?(?)*-*  r(f)c^ 


where  in  the  last  equation  it  is  clear  that  (0„   is  the 


N 


surface  of  the  N- dimensional  unit  sphere  and  P  is  the 
well  known  gamma  function.   On  the  other  hand 
for  N  odd  we  write  using  2.23 


2,25 


To  evaluate  the  integral  on  the  right  of  the  last  equation 
it  is  necessary  to  make  the  following  change  of  variables* 
Let 


2„26       e  =  z(f)^cv^))"' 


13 


then 


2.27  ^  =  ^f)(v^° 


and 
2.28 


i  _l  r  _  , 


Hf)lZ 


A^-^0         0  I 


1     r00    -J 


-2cf)"*  r(i). 


The  value  of  the  above  Integral  is  obtained  by  direct 
calculation.   Substitution  of  2.28  into  2.25  now  yields 

2-29       fe  il^f^--i^(.^m^)..ir(i)  Q„ 


=  inp)(n)i   p(^)g 


where  use  has  been  made  of  the  well-known  recurrence 
relation  for  the  gamma  function.   The  desired  result  is 
now  obtained  by  noting  from  page  223  of  [2 J 


2.30  0N-ZTTK[rCi)} 


Ik 

and  substitution  of  this  expression  into  2.2l|.  and  2.29  yields 

2.31       i^  ff  V***  Js  =  <P  (*.>*,} 

for  N  either  even  or  od 

3.   An  Auxiliary  Fun      for  the  Parabolic  Problem 

Here  our  aim  is  to  construct  a  sufficiently  smooth 
function  r   which  will  be  employed  in  the  following  chapters 
to  obtain  solution  bounds  for  the  parabolic  problem.   The 
desired  function  -j1   will  be  defined  so  that 

2.32  -<^T  <  oo  for  (x,y)6V 
and 

2.33  J{^)<0  for(x,y)£V 

where  J  satisfies  the  conditions  of  the  previous  section. 
In  addition  to  the  above,  it  is  necessary  that  "Z7  have  a 
singularity  at  peD(y_),  which  is  of  a  slightly  lower  order 


o 


than  the  singularity  of  the  parametrix  )&   p  at  p£D(y  ). 

The  desired  function  is  constructed  by  first  considering  for 

N  £1 

and  computing  for  (x,y) £  V 

2.35  W)"(s.-«jr1*J  [-4  \^->)ict)a  tfteVXa.-**) 


15 


H-H 


U^)2(<jo)o\^)(^ 


*& 


^* 


+(?)*(Kiy-hteN-)**i*Y  *  %t%-&y 


-  Hf  [-  few- ')  *  (f>)  f>  \8K%-  *£) 


i    ^ 


h  ^ 


-^/,)"\2Ww)^(t)>^(w)2(•^vw)^^(1l)^^CxV)%jp)(><*-o%cfJ^1^ 

o 

with  g1J,  gij,  h  and  h  defined  by  2.19  and  2.21. 
For  (x,y)6V  and  sufficiently  close  to  p,  the  second  term 
in  the  brackets  on  the  right  of  the  last  equation  dominates 
the  first  and  third  through  seventh  terms  while  the  eighth 
dominates  the  remaining  terms.   Due  to  the  boundedness  of 
the  g's  and  h's  and  the  piecewise  continuity  of  a^  in  V, 
it  follows  that  there  exist  positive  numbers  g^  and  (Dq   such 
that 


16 


2.36 


rj&)u 


j(qr)<o 


Let  o^_  be  such  that 


2.37 


^* 


%^-«f>  W)(*-D*  -  <*»-)*<#*   e  <o 


for  all  (x,y)e|x,y  (x,y}£V,  o s  yQ-y <  S^}  .   We  set 


2.38 


S0   =  min  {  8  ~  ,  S0l] 


and  then  define  the  desired  auxiliary  function  as  follows  in 
terms  of  non-negative  constants  <*  and  tr 


r 


*  +oc  H^'i)     ^t'  (tv-i)  scp)/?  (?«**-$'} 


-<Z 


(x3^)€V       3      00^-3    <  St 


•39     ^r J  =  si * '""**  (^^X^l-^^-o^K^^^)]"'} 


4 


17 

The  constants  cC  and  xr    appearing  In  the  last  equation  shall 

be  explicitly  chosen  in  the  proof  of  the  theorem  which 

follows  this  paragraph.   On  the  surface  D(y0J-{pJ,,  the 

function  -j7   and  its  derivatives  are  defined  as  the  limits 

of  these  functions  as  (x,y)  approaches  (xs y  )  with  (x,y)  £  V 

o 

and  x-#x  ;  thus  we  have 
'  o 

2.^0  ^U,y0)  -  hm    ^Fu9j)  x^xQ 


y+y 


0 


2.1+1         %L^»yo)  =  ^  %(*SJ)  x*x0 


o 


2.1+2         «f,,s(x,y J  =-2im¥^  .(x,y)       x* 


y*y0 


X  1=    X, 


and 


2. [j.3         -^—       =  JLixa.  fLL 


^y      y-»yrt  "^y 


x^xf 


•o 


Theorem  2.2,   The  function  ^?  (o<.  and  fcr  appropriately 
fixed)  is  twice  piecewise  continuously  dif ferentiable  in  x 
and  piecewise  continuously  diff erentiable  In  y  for  (x,y)£V 
except  at  p  =  ^o'^o^  an(^  ^s  SUC^L  that 

2.1^         ^O  ,  !(¥)&o  (x,y)€V 


and 
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2.1^        czf>o   ,  J(T)<o    (x,y)£{xsy|(x,y)€  V,o<  £<yo-y} 

where  the  coefficients  of  J  satisfy  the  requirements  of 
section  2  of  this  chapter  and,  in  addition,  the  condition 


2.I4.6         min  Z  -  mp  >o 
V 


Proof:   The  function  ^V   obviously  satisfies  the  continuity 
and  differentiability  requirements  of  the  theorem. 

Our  aim  shall  now  be  to  establish  2-»\\\\   and  2.d\.%   for 
properly  chosen  constants  ^C  and  at  .   Prom  2.36,  2„37  and 
the  definition  of   r  it  is  clear  that  the  desired  in- 
equalities are  satisfied  for  all  o(£   o  where  o  ^  y0~y  £  §0 
and  P^ft  0°   For  P^     and  °  ~  J  "*y  ^  S  we  compute  by 
means  of  2.39 


z.ki  im  -  t c*)  + « (a^y  ^J  { §| (^ 


+  (^)<^-tf)?  ~  b»-i)w$/£W  *} 


•  ^f(^-0^^SA/(^))"'j  . 
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For    (x,y)  €-V   and  y  »y  >S       we  have 

o     o 


248  jcr)  .  ^'^i-toft*  (M^u^o  C*'V-W" 


-  (vw)-'ft'>(2M-o «(p)a^(f)  t^y" 


+*  1 If  Va-)*" "  *  l  (/3o-  i)1 " 

Having  calculated  the  above,  we  are  now  in  a  position  to 
choose  the  non-negative  numbers  oC  and  J(r .      Since  2„37 
is  valid  in  the  domain  of  2„!|7.*  it  is  clear  that  o£  >  o 
may  be  fixed  so  that 

2^9         J(*f  )£o       (x,y)6  {xsy|  ^0,oiyo-y  *  £  q] 
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and 


2.50  J(^)  <  o  (x,y)e{  x,y|  p>f>Q,o<  Z  <  yQ-y  <  S0] 

By  Inspecting  2.1j.8  it   is   clear   that    ir  ^  o  may  be   chosen  so 
that     • 

2.^1  !&)<o  (x,y)^{x,y  |(x,y)6  V,   y  -y  >  S    | 

*-  o  OJ  ' 

In  order  to  establish  that  the  auxiliary  function 
is  non-negative  in  V  and  is  strictly  positive  in  every 
closed  set  of  V,  we  need  only  consider  equation  2.39. 

Subject  to  the  foregoing,  it  is  clear  that  the 
theorem  is  established;  however,  as  a  final  remark,  we 
note  that  for  certain  differential  operators  we  may  choose 

2.52  %    ^  (x,y)e{x,y|(x,y)£V,o<yo-y  <  Sj 

since  J(4")  satisfies  2.I4.9  and  2.£0.   When  this  condition 
is  applicable,  one  would  of  course  choose  oC  =  o  in  2.39. 
The  conclusions  of  Theorem  2.2  shall  be  employed 
in  the  next  chapter  to  obtain  pointwise  bounds  for  the 
solution  of  non-linear  normally  parabolic  problems.   In 
the  first  section  of  Chapter  IV  the  function  3-  =  Z(x,y) 
is  not  assumed  to  be  strictly  positive  and  2. [(.6  is  replaced 

by 
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2.53         min  a  >  o 

V 


and 


2. ft         Z(x,y)  >o  (x,y)£  D(yQ) 


Under  such  an  assumption  it  is  obvious  that  Theorem  2*2 
is  not  valid  since  we  would  in  general  be  unable  to 
choose  9r   so  that  2.[j.8  defines  a  strictly  negative  number 
Because  of  this  reason  we  use  the  fact  that  a1^  is  a 
positive  definite  matrix  and  introduce  a  slight  different 
auxiliary  function.   This  new  function  is  much  more 
difficult  to  construct  since,  besides  having  a  more  com- 
plicated form,  it  is  necessary  to  calculate  a. . (x  ,y) 
for  o<y-$y0.   Thus,  even  though  the  revised  auxiliary 
function  could  be  used  in  the  linear  normally  parabolic 
case,  the  function  as  given  by  2.39  would  be  easier  to 
construct  since  the  inversion  of  the  matrix  a1^  is  only 
necessary  at  the  point  p  £  D(yQ) • 

When  conditions  2.53  and  2„5i|.  replace  2.!|6,  we 
define  h*    as  follows: 

2-^  ^  =  (3.-  *)"  ^  *?  H»-d  W  V^-*)P 
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for  (x,y)£  V  where 

2.56  p   -  *± ^  (x0,y)  (x1-**)  (x  j-x^) 

For  our  present  definition  of  ^  we  notice  that  2.35  li 
valid  with 

2.57  /O   replacing  xo 
and 


2.58    a.      (x  ,y)     replacing    0,  (p) 
rs  o  rs 


In  the  revised  form  of  2.3£  it  is  also  necessary  to  add 
the  following  term  on  the  right  side  of  the  equation 

2.59       4**l^few-.)(8«(V^)'*1)" 


By  exactly  the  same  arguments  as  previously  used,  we  deduce 
that  there  exist  positive  numbers  2  £  q-^  and  Aw  such 
that 


and 
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2.60  J(<n  s<o  (x>7)e{x,Y\/o$fQ1,oi7o-7S  2£Q1] 


2.61  J(^)<o  (x,yK[x,y    ^'£/}01,oiy0-yS25:o:Ji 


Since  Z(x,y)  is  strictly  positive  on  the  surface 

D(y  )  it  follows  from  the  conditions  of  the  first  section 
17  o 

of  this  chapter  that  there  exists  a  positive  number  £ 
such  that 


02 


2.62 


l^v-1^ 


£  <  o 


for   all    (x,y)f{x,y  |(x,y)  6V,o<yQ-y<  25Q2j 


We    set 


2.63 


So  =  rain{S01,  S02] 


and  notice  that 


2. flj.         J(5  )  «  o       (x,y)€  [x,y|(x,y)£V,o^y0-y<2<S0] 


and 


2k 


2.65  J(oC)<  o  U,y)e{x,y  |(x,y)6V,o<£  i  y0-y£2SQj 


where 


2 .  66  <*  =  oc  (^y    *  a^o  {-ft*-,)  *(^KEA/^o-^))"'] 

Since  2.62   is   true,    it   is   evident   that   2.6fy   and   2.65  are 
valid  for  all   dC>  o.      By  comparing  "^T  and  o(   ,    we   observe 
that  we  may  choose  ot   -  °<-  ^,    such   that 


2.6?  JC5^  )  So  (x,y)£{x,y  |(x,y}*  V,o$y0-y$2§0] 


and 


2.68  JC=Ffo(  )  <  o  (x,y)e{x,y  |(x,y)  €■  V,o<ti  y0-ys2S0} 


Next  we  introduce  the  function  ^<  given  by 


or 


2-69         $--  is1  -cv^)"}- 


•  a^jd  {-  (sws.y'/tf  O-/)  *(f3  /3"\? 


+-££■  to? 
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where  for  the  present  /3  is  an  arbitrary  non-negative 
number.   To  fix  this  arbitrary  number  we  compute 

2.70  x(e)  =   {  s;1  -c^-^y}- 


For  &  <  y  -y  £  2  §  and  (x,y)  £V?  the  last  term  on  the 
o   o      o 

right  of  the  above  equation  is  a  strictly  negative 
number.   Next  we  consider  the  following  limit 
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2.7i  j&^       j(p)    --    {S;'-C^y] 

•f(ys0y'^  (M»-0  ts-t/o)   +  If] 


By  the  remark  following  2.70  and  due  to  the  form  of  2.71, 
we  deduce  that  there  exist  positive  numbers  /^np  and  ($ 
such  that  for 

2.72  P    =fii 
we  have 

2.73  J(#)  <o       (x,y)e{x,y|^^02,o^yo-So}( 

We  are  now  in  a  position  to  define  the  desired 
auxiliary  function.   As  was  mentioned  above,  this  function 
will  be  used  when  bounds  for  the  solution  of  a  degenerate 
problem  are  sought.   The  function  is  explicitly  given  by 
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=<¥  + 


--<% 


ct 


+  (*W  sjT  &*-j)  *(fi)/3,  /°.i  +  '  j ' 


Z. 

or 


4^  +cL 


f^oi^^^-K 


On  the  surface  D(y0)-jpj  ,  the  function  -/?  and  its 
derivatives  are  defined  as  the  limits  of  these  functions 
as  (x,y)  approaches  (x,y  )  with  (x,y)6V  and  x^x  .   We 
now  show  that  with  <C  properly  chosen  the  function  ^f7" 
satisfies  the  conclusions  of  the  following  theorem. 

Theorem  2.3.   The  function(^~,  defined  by  2.7^, 
is  twice  piecewise  continuously  diff erentiable  in  x  and 
piecewise  continuously  diff  erentiable  in  y  for  (x,y)6V 
except  at  p  =  ^o'^o^  and  ^s  SVLC^i   tha'fc 
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2.7S        J(<2F)^o  (x,y)&V 


2.76         J(<sf)<  o       (x,y)e{x,y|  (x,y)€  V,o<£*  yQ-y} 


and 


2.77         maxj-^j 


M 
V 


where  M  is  a  finite  number.   The  coefficients  of  the 

differential  operator  J  satisfy  the  requirements  of 

section  2  and  in  addition  the  condition  that  the  partial 

derivative  of  Ou      (x  ,y)  with  respect  to  y  is  piecewise 

rs  o 

continuous  on  the  line  (xQ,y)  £  V.   The  function  S  is  such 
that 


2.78         min  %  ^  o 

V 


and 


2.79        Z(x,y)>o  (x,y)6D(yo) 


Proof:   The  continuity  \nd  differentiability  of  ™  are 
obvious  for  (x,y)e{  x,yj(x,y)GV,  o£yQ-y£  S  J  .  For  the 
lower  part  of  V  we  calculate 
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2.80 


/°  "V  0  2 


«tf  fft'uS»)7J.^w-,)*^)/j»  }  +^  +<* 


and 

2.81 

'H^j'^  (iM-,)a(f)  *./-o^)  (**->>*)}' 

=  Aw.     i^J. 

Our  task  is  next  to  select  the  number  JL   so  that  "/^ 
satisfies  2.75  and  2.76.   Prom  the  remarks  associated 
with  equations  2.5>5>  to  2.73  it  is  clear  that  these  in- 
equalities are  satisfied  on  the  set 


30 

For   Z)  -^02   and   0^3r^yo~  &o  we   comPute 

2.83  Tfr>  { v-c^rl  {(*«-'  av-o*(^K-  A^^f 


'**t»  {"  (fiu&fl,  t2»-}*tp)fc}  +  C-  v$ 


-i+ 


•fe1-  ^i  ifH-^^rvo*^^)/? 


+ 


2. 


+(*»sjy  (z»->)  «<>)  ^  /*.»  +  /  i 


,J^fi-(^£0)"'/5,  (iw-j)*(^3  +  ?(?+*). 


In  order  to  determine  an  appropriate  choice  of  the 
constant  tL   we  note  that 

2.81+        ai0'airC*o,y)  W-aJPSj  lVy)  (x»-x«)£  aajl/)  ^ 

for  all  (x,y)  belonging  to  the  domain  of  definition  of 
2.83;  hence  it  is  evident  that  0L   may  be  chosen  so  that 
2.83  defines  a  strictly  negative  number. 

With  d.   so  chosen  it  is  obvious  that  the  finite 
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number  M  given  by 

2.8£         M  =  maxf  -qF] 

exists  since  T7  is  bounded  everywhere  on  V  except  at  the 
point  p  and  at  this  point  one  has 

2 .  86         JU^u    ^T  >,   o  , 

In  future  investigations  bearing  on  the  subject 
matter  of  this  thesis,  an  attempt  will  be  made  to  replace 
2.79  by 

2.87  B(p)  >o. 

Next  we  shall  establish  the  existence  of  an 
integral  which  will  often  appear  in  the  bounds  which  are 
developed  in  the  following  chapters. 

Theorem  2.1|.   The  function  ^^  is  such  that  the 
integral 

2.88  -US[*eF)]"[3(.XftfJlv 
exists . 
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Proof:   As  in  the  other  theorems  of  this  section,  the  proof 
is  completed  by  discussing  two  cases.   When  2.J|6  is  satisfied, 
it  follows  from  2.35  that  there  exist  positive  numbers  S  -, 

and  /O.  which  are  such  that 

N+3 

2.89   -J(£)>  (y  -y)  2 


|z(p)(y0-y)  +  (32N2)"1(2N-l)Z(p)2^' 
•exp{  [-  SNCy^y)]-1  (2H-l)Z(p) ^21\ 

for  /o  */°-\i  ft      and  oi  Y  ~7i  S-,1  5  •   With  ^o  2  replaced  by 
/£>   it  is  evident  that  a  similar  result  follows  for  ^A  in 


case  2.79  is  applicable.   Hence  we  have 

N-fl 
2.90  -Iffi)}  |(y0-y)""^"z(p)exp{-  [8N(yo-y) 


-1 


(2N-l)Z(p)^2J} 
for    P  9ft    i  ft  .   <  ft      and  o$y-y<    £    <    S    .      Next  we   consider 

2-91         Jff[-?fr)V[T(*JlJ.v 

v  'J 


y-y 


v/ 


where 


2.92  V1  =   {(x,y)      /05/o    ^/?1,    olyQ-y<    o^}. 

The  first  integral  on  the  right  of  2.91  clearly  exists  since 
its  integrand  is  bounded  over  the  domain  of  integration; 
also,  since 


2.93 


(8N)"1(2N-1)  <|  , 
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the  second  integral  on  the  right  of  the  aforesaid  equation 
is  seen  to  be  finite  by  application  of  2*90  and  Theorem  2.1. 

3.   Generalized  Form  of  an  Identity  of  Hbrmander 

The  identity  under  consideration  was  derived  by 
Hormander  |l6J  for  the  purpose  of  obtaining  a  uniqueness 
theorem  and  estimates  for  normally  hyperbolic  partial 
differential  equations.   Essentially  the  same  identity  was 
employed  by  Payne  and  Weinberger  |_2lJ  to  obtain  solution 
bounds  for  second  order  elliptic  equations .   In  what  follows 
a  slightly  generalized  form  of  the  same  identity  is  derived. 

Let  f1(x,y),  .  .  .  fN(x,y)  be  any  set  of  functions 
which  are  piecewise  continuously  diff erentiable  in  x  for 
(x,y)6  V.   By  direct  differentiation  and  application  of  the 
divergence  theorem,  we  have  for  any  function  ^  which  is 
continuously  diff erentiable  in  x  for  (x,y)  G  V  and  whose 
second  derivative  in  x  is  continuous  in  the  interior  of  a 
finite  number  of  subregions,  the  sum  of  which  is  D(y),  where 
y  is  such  that  o  <y<y   (in  the  above  we  suppose  that  D(y) 
is  a  domain  of  N  dimensional  space) . 


a 


3k 

The  functions  a1  J  ■  aij*(x,y)s  i,j  =  1,  2  .  .  .  N,  which 
appear  in  the  last  equation,  are  assumed  to  be  symmetric, 
piecewise  continuously  diff erentiable  for  (x,y)6V  and  such 
that  2.7  is  satisfied.   (For  the  results  of  this  section  the 
piecewise  continuity  of  the  derivative  of  a  ^  with  respect 
to  y  is  not  needed^  buty  since  this  condition  shall  later  be 
necessary,  it  is  included  here.) 
Let 

2.95       Tb~  ^  %i^^  (*»«&)€  5 

and  then  note  that  the  vector 

a.96       (^H*%rVU>V/b) 

defined  almost  everywhere  on  S    is  orthogonal  to  (n.,  n  ) 

where  in  2 .96  we  have  made  use  of  2.5  and  2.7 •   Prom  the 

above  it  is  clear  that  (7*^»°)  *s  a  tangent  vector  to  the 

surf ace S  ,  is  perpendicular  to  the  vector  (0,1)   and  is 

defined  almost  everywhere  on  5  .   Following  closely  the 

procedure  of  section  2  of  I 21  J  we  renormalize  by  setting 

1  JL 

2.97  (^,0)  «  crH^Vq]2^^^8]  2>  0) 

where  a.  .  is  the  inverse  matrix  of  a1**.  For  (x,y)6<5 
we  have 

2„98        a^t1^  =  a  Jn.n.  =  n 

where  n  is  defined  by  this  equation.   Due  to  the  above,  the 
directional  derivative  on  S, 
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2.99 


is  in  a  tangential  direction  which  is  perpendicular  to  the 
ny"  direction..   For  (x,y)6  5  we  notice  that 

2.100 


^v.vw«^ 

=r* 

s^* 

^        ■ 

\ 


Also 
2.101 


^'^-^ 


2  /»v//>y 


>  ^ 


^  ,  /  /*JT)*- 


9      * 


Hence  almost  everywhere  for  (x,y) £  5 


2.102 


4y  „       -/r^V^^i^/^  fcl    -'[miV^f^ 


■krrpJ: 


"A  A' 


-/ 


(#(H/Vr? >9 


/n 


-/ 


5v/+VW 


Moreover  almost  everywhere  on  5  we  have 
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2-103  ^A'l'vV'f/' 


V'*** 


1^  ^*J 


=/»,fY,  - 


^w 


V   7 


tf7~"t- 


J/i> 


^ 


and   thus 

2,101+  f  «».  3V    +  fyJt'St  ^  nf    ?>-  • 

Using  2.102  and  2ol0i|  we  rewrite  2.9l|.  as 

rrrt*  A*  ^  **  />*  *'*  ?k  -^1 ,       < 


We  now  choose  f  ,  k  -  1,  2  .  „  .  N,  so  that  a  bound  for  the 
square  of  the  integral  of  the  conormal  type  derivative  2.95 
may  be  obtained  from  2«105«   Previously  we  required  that 
these  auxiliary  functions  be  piecewise  continuously  differ- 
entiable  in  xj  we  now  further  assume  that  r  n, ,  is  bounded 
and  has  a  positive  minimum  on  5.   Such  a  set  of  functions 
was  considered  by  Payne  and  Weinberger  in  section  3  of     21 I. 
Prom  this  reference  we  note  that,  if  $  is  starshaped  in  x 
with  respect  to  some  line  perpendicular  to  the  x  plane,  then 
we  can  choose  this  line  as  the  x  origin  and  take  f  =  x  . 
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Also  if  V  is  a  domain  between  two  boundaries  each  of  which 

is  starshaped  in  x  with  respect  to  the  above  perpendicular 

line  and  if  V  contains  a  cylindrical  shell  of  radius  r_ 

o 

Xc  k 

from  this  line  then  we  may  take  f  =  x   (r-r  ) »   For  further 

information  bearing  on  this  choice  of  the  auxiliary  functions, 


the  reader  is  referred  to  Appendix  A  of 
defined  above,  set 


[l7 


With  fx  as 


2.106 


m0  =  min  n 
3 


-1 


"Ky. 


Since  a*  3  is  positive  definite,  there  exists  a  constant 
Ci  which  satisfies 

2.107  -  [f^aij-f,yk-f,yk  +  ***.JJ]w«;  x\  rti 


for  all  (x,y)  ev.  (C,  is  an  upper  bound  for  the  largest 
eigenvalue  C,  of  the  coefficient  matrix  on  the  left  with 
respect  to  a1J0)   For  a  crude  upper  bound  for  C-j_  we  define 


2.108 


Jr 


iJ  =  - 


f,Vj-^Sjk-^Jaik  +  fka>ij 

.  k      k      k         k 


then 


*See  also,  Bramble,  J.  and  Hubbard,  B.,  Some  Higher 
Order  Integral  Identities  with  Application  to  Bounding 
Techniques,  in  print . 
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2-109        7_  ,   khjll  -  £%!*_ 


*x? 


where 


2.110        B  =  maxOj^)' 

V   *j 

Substitution  of  2.106  and  2.107  into  2.105  yields  after 
applying  the  indicated  inequality  to  the  resulting  expression 


2.111 


where  c(  is  an  arbitrary  positive  number.   By  setting 
2.  112       d  =  T 
we  may  write 
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where 


V    J* 


2.114       %^3{K^\)'\srt^tfM 


and  clearly 

2.115        ^  -  0 

If  V  =  o  for  (x,y)6S  -  In  the  following  chapters  great  use 
shall  be  made  of  2.113  when  the  boundary  value  problem  under 
consideration  is  of  the  Dirichlet  type. 


CHAPTER  III 
NON-LINEAR  NORMALLY  PARABOLIC  PROBLEMS 

lo   The  Non-Linear  Mixed  Problem 

In  this  section,  bounds  for  the  solution  function  of 
a 'boundary  value  problem,  which  is  normally  parabolic,  will 
be  derived  where  the  solution  function  is  initially  given 
and  thereafter  satisfies  a  non-linear  Neumann  condition  on 
the  boundary  and  a  given  differential  equation  in  the 
interior. 

We  assume  that  V  is  an  N  +  1  dimensional  domain  with 

boundary  D(y  )U  SUD(o)  which  is  such  that  the  conditions  of 
o 

the  first  section  of  Chapter  II  are  satisfied.   Our  goal  is 
to  determine  pointwise  bounds  for  the  solution  of 

J(w)  =  (<l*wJ^-  *  V$  =  f  (*^  w>  v w)  b*$ev 

W(>30)  =  &(*)  (*4)6D(*) 


3.1       4 


—  oO  <    W  <  OO 

where  W  is  continuous  in  VUD(yQ),  continuously  differentiable 
in  x  for  (x,y)  6  V  and  is  such  that  its  second  derivative  In 
x  is  continuous  in  the  interior  of  a  finite  number  of 
subregions,  the  sum  of  which  Is  D(^),  where  y  satisfies 

ko 


in 

o  <y<yQ»   Also  we  suppose  that  the  first  derivative  of  w  in 

y  is  continuous  in  the  interior  of  a  finite  number  of 

subintervalss  the  sum  of  which  is  the  line  (xsy)  for  every 

fixed  x  where  (x,y)€.Vo   In  addition  it  is  necessary  that  the 

solution  of  3«1  satisfy 
r 


3-2         < 


lim  w(x,y)  =  g(x) 
(xsy)->(x,o) 
(xsy) a 

lim  a*J(x,y)w8  (x,y)  =  ai^(x,y)ws  (xsy) 
(x,y)-^(xsy)     1  ,     x 

(x,y)  €V  -*,=  /3£3 A/, 

(x,y)£S  ^ 

The  boundary  value  problem  is  assumed  to  be  such  that 
the  components  of  the  symmetric  matrix  a1^  =  a*J(x,y)  are 
piecewise  continuously  differentiable  in  x  and  are  such 
that  2.7  is  satisfied  in  V.   In  the  closed  set  V  the  function 
Z  =  Z(x,y)  is  piecewise  continuously  differentiable  in  y  and 
satisfies 

3«3         min  Z(x,y)  =m  >o. 

The  boundary  data  T   ,  JL    and  £  are  assumed  to  be  integrable 
and  square  integrable  over  their  domains  of  definition  for 
any  bounded  continuous  function  w  which  is  continuously 
differentiable  in  x  for  (x,y)  £V.   In  addition  f  and  JL 
are  required  to  satisfy  a  Lipschitz  condition  in  all  but 
their  first  two  arguments.   Hence  there  exists  positive 
numbers  M- ,  MpS  M  ,  i  =  1,  2  .  .  .  N  such  that 
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34 


IfO 


s 


l-£fr,*3  vO  -  ^  (x,  ^  wO  *  f\  |  W,  -  Wz  | 


(x,^)€5  j,  -«o  <w<   oo  # 


Except  for  the  fact  that  Jt  =  £(x,y,  vy)  is  required 
to  satisfy  a  Lipschitz  Condition  in  w,  the  problem  of  this 
section  is  a  generalization  of  a  problem  whicn  Friedman 
considered  in  [ill .   In  that  paper  the  elliptic  part  of  the 
differential  equation  was  the  usual  Laplace  operator,  T 
was  identically  zero  and  Ji   was  a  monotone  non- increasing 
function  of  w  .  Under  these  restrictions  the  solution  of 
the  boundary  value  problem  was  shown  to  exist  by  means  of 
the  Schauder  Fixed  Point  Theorem  [26J.   In  a  concluding 
remark  the  author  indicated  that  in  a  later  paper  he  would 
attempt  to  extend  the  existence  proof  to  general  second 
order  parabolic  equations.   Whether  or  not  such  an  extension 
has  been  carried  out  is  not  known.   The  existence  of  a 
solution  for  a  similar  problem  was  however  established  by 
Piskorek  [22]  and  also  by  Pogorzelski  [23] • 

Let  p  GD(y  )  be  the  point  at  which  a  bound  for  the 
solution  is  desired.   We  assume  that  the  functions  a1 J  and 
2  satisfy  a  Lipschitz  condition  at  p€D(y  )  for  (x,y)  £  V, 


hi 


that  is,  we  assume  that  2.8  and  2.9  ai*e  satisfied.   In 
addition  it  is  necessary  to  suppose  that  f(xsy,wsvw)  is 
bounded  in  some  neighborhood  of  the  point  p  which  is  con- 
tained in  V  for  any  continuous  function  w  which  is  continu- 
ously diff erentiable  in  x  for  (x,y)  6  V.   Hence  for  any  such 
function  w  there  exists  numbers  8  ,  p   9    and  & (w)  all 
greater  than  zero  such  that 


3»5 


f (x,y,w,  v  w)  ^  %-(w) 


for  f>  if>2   and  y  -  t>2$jiy   . 

The  desired  bound  for  the  solution  is  obtained  by 
choosing  an  arbitrary  function  lp  (x,y)  which  is  such  that 

K  ID 

JC1^)*  5v  and  l^  (x,o)  approximate  f(x,y,  ^,^^),  ^(x,y,  ^) 
and  *  (x)  respectively  where  we  further  require  that  (P  be 
twice  piecewise  continuously  differentiable  in  x  and 
piecewise  continuously  differentiable  in  y  for  (x,y)  £  V. 
Let 


3.6 


^U,y)  -  w(x,y)  -  <P(x,y) 


and  compute 


3.7 


where 


< 


(^^S  J  -cOitf^toQ 


hk 


3.8  F(x,y)    -  f  (x,y,  ip,  7tf)    -   J(  </>) 

and 

n  ^<P 

3-9  L(x,y)    =  /(x,y,  4> )    -  ^    . 

Due    to   3' 3  an<3   3«k  we  ma7  write 

3.10  JCmSjIy/   +-M^J    +    F(*/*) 

and 

3.11  r£  <  f\)<fj  +  L  (*,/*). 


With  7(f)    and    V       defined  by  2.6   and  2.10   respectively, 
from  the   divergence   theorem  it  follows    that 


•s(a)  *  r  we)  f 


where 


3-13  V(y)    =  j(x,/w)|  (x,/>j)  £D(^  ),    o</»<yJ 


and 


3.14  5(y)    =  {(x,^  )|  (x,/n  )£  D(/m  )-D(/H),    o<^<y] 
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In  3* l4  we  let.  y  approach  y  and  use  Theorem  2.1  to  obtain 


o 


which  becomes  after  substitution  of  equations  3*7  to  3°11 

3.16  lw(*o^)  -  7(f>)l 

+{ f /is}4  {<{f  (#  *y  tfi^+  i$ry«^} 

where 

3.17  X(f)=-  if(*orf„)  +   ttifr^Vpds 

bio)     r  *-»*o       V(g)       7 

The  terms  on  the  right  of  3' 17  are  all  known; 
however,  it  is  necessary  to  prove  that  the  last  integral 
converges .   For  this  purpose  we  write 


k6 


=  i^w    {  ///   V.FJv  f  J//    X„F<lv] 

n>  muss,)     vmv$£) 

v-v(ft5sj  v-v^A) 


where 


3-19  v(# ,  O  =  { <?,$  \f>  W  ,  o  -<  ^  <  St  } 

and  on  the  surface  Dfy^-jpl  it  is  clear  that  tf   is  defined 

as  indicated  in  Section  2  of  Chapter  II.   We  now  consider 

the  three  integrals  on  the  extreme  right  of  3-l8»   The 

first  exists  since  o        is  bounded  in  the  domain  of  inte- 

P 

gration  and  f  is  inte/^rable  over  V.   The  second  integral 
exists  since  its  integrand  is  bounded  over  the  domain  of 
integration.   For  the  last  term  we  use  the  fact  that 
[f  (x,y,  if  s  9$  )-J(  <f  )]   is  bounded  over  V(^2,  So)      due  to 
3.5  and  the  definition  of  J;  thus  the  integral  exists  since 
^f      is  easily  seen  to  be  positive  and  integrable  over 

v(/*2,  82). 

Prom  the  above  remarks  it  is  clear  that  all  the  terms 
on  the  left  of  3.l6  except  w(xQ,y0)  are  known  and  hence  to 
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obtain  a  bound  for  the  solution  it  is  only  necessary  to 
estimate  the  expressions  on  the  right.   To  obtain  this 
result,  we  introduce  the  function  T  defined  by 

3.20  f  -  £exp[-K(y0-y)J 
and  calculate 

3.21  JCM  =  J(?exp[-K(yo-y)])  =  J(£)  exp  [-K(yQ-y)] 
where 

3.22  J(?  )  =  J(  f)-Af  ? 

Let  f  ,  i  =  1,  2  .  .  .  N,  be  the  set  of  auxiliary  functions 
which  was  introduced  in  Section  l±   of  Chapter  II,  that  is, 
the  f  are  piecewise  continuously  differentiable  in  x  for 
(x,y)-6  V  and  are  such  that 


s 
With  the  f*  as  thus  defined  we  may  write 

~£  /    -Iff   ~^ 
5 


3.23         o  <mi  =  min  f  n^. 


3.24  S/vUs  *^tf<fzfA*zJ.s 


^usujUv  +±sff?fU* 


<<n< 


SfU%^^^)Sff?% 


where     °C    Is   any  positive  number  and 
3.25  Mx  ■  maxfZfVj 


^8 
3.26  M~   =  max  \tA  7 

In  arriving  at   3°2i^  we  have  made  use   of 

3-2?         (f\f   '  ^^  ^*4 
which  is  easily  established  by  expanding 

3.28         f  cr?,  -{**>,  f. 

From  the  divergence  theorem  it  is  clear  that 

3.29  US*H*)JL*a-SSS*u**Hih  +11"*  ^ 

Using  3»7  and  3°  20  to  write 

3-30  ft*)*  ^  lH+Pll\%J  +  f^)«f[K(3,l)] 

and 

we  have  by  means  of  the  above 
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3-32         tisf  9  iwjivU  summtUv 


and 


3-33 


f  ^  XT  L4^DpLK(^-^)J  ,/5 


where  ct  ,  oC   ,    and  °^i  are  arbitrary  positive  numbers  and 
3-34 


/V 


M,  =Z(M1)2, 
3    / 


Substitution  of  3.32  and  3.33  into  3-29  yields  after  the 
necessary  transposition 


where 

3.36 

and 


i<3  w  f%f)N*-a)]Jv+  wvf ^k-^5"^ 


M.     =  max  L  t)Ax  j 
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3.37        M^  =  max  {-n  zj. 
->    s     * 

Except  for  the  first  and  last  integrals,  the  right  side  of 

3.35  is  completely  known;  the  last  term  may,  however,  be 

neglected  since  Z  satisfies  3»3'   In  addition  we  note  that 

the  left  side  of  3-35  has  the  same  form  as  the  right  side  of 

3.2I4.;  and,  after  a  close  scrutinization  of  these  inequalities, 

it  is  evident  that  they  yield  a  bound  for  the  integral  of  the 

square  of  f    over  S  provided  the  arbitrary  positive  constants 

oCj  and  K  are  properly  chosen.   In  a  particular  problem  the 

best  choice  of  these  arbitrary  constants  would  depend  on  the 

functions  P,  G,  and  L  and  the  M's  appearing  above.  For 

definiteness  we  assume  M  ^o  and  make  the  following  selections. 

(If  Mo  ■  o,  obvious  modifications  are  made  to  obtain  the  same 


result) . 

Let 

3.38 

% 

a 
0 

"M3 

3.39 

h 

=      1 

3-1*0 

1 

aomlf 

34l 

*3 

=  1 

-1 


and 

1     J\  M3  1 

342         *-=5vf  **  +  4;  +  5 


'2(2M2+l+Mc;)M1         M^ 
+  (2M0+1+Me)  { ?-2 — =  +  —  >  V    =  K 
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where  we  have  assumed  that 


343         %  +  |  +  -|  M^>o. 


If  this  is  not  the  case,  then  one  need  only  choose  °C  i  so 
that 

^    °C)i    i 
3.kk  M2  +  2  +  I  M5>0 

and  the  same  result  is  now  obtained  by  appropriately  defining 
K  and  cC-^.   Substitution  of  the  above  into  3-35  and  3.21^ 
gives 

34*        SSVUs  *  ktii+l+M${SSG***f(U,>j0)Js 


Ho) 


v 


tff  f2  >^  [z  k,  Oj,-$}j  v*{h\[z  h,  fy-*)]*} 


where  K-,  is  defined  by  34^  (assuming  343  is  valid).   We 
use  the  last  inequality  and  3 °20  to  obtain  if  K^ ^  0 


+ 

V 


and  if  %.,%  0   then 


JXf  flMf(-ZK,^Jy  4  Ml* tff&l  %)*!*} 


347  #Ws  «Hm4li|f^(«^ 
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The  right  side  of  the  applicable  inequality  is  known  since 
Ps  L  and  G-  are  square  integrable  over  their  domains  of 
definition;  hence  we  may  write 

where  p  is  a  known  number  depending  on  the  function  $   . 
After  having  calculated  3<>4^>  we  note  that  a  bound  has  been 
obtained  for  the  last  term  on  the  right  of  3»l6. 

If  a  function  I  p  can  be  determined  which  is  such 
that  the  conclusions  of  Theorem  2.1  are  valid  and  in  addition 
the  condition  that 

3.I49       J(fp)  -  0, 

then  it  is  clear  that  lp  is  a  fundamental  solution  of  J  and 
the  first  term  on  the  right  of  3«l6  is  identically  zero. 
Dressel  [3]  and  [lu  established  the  existence  of  such  a 
solution  when  the  coefficients  of  the  differential  operator 
are  sufficiently  smooth.   Prom  these  references  it  is  clear 
that  even  when  its  existence  is  guaranteed  such  a  function 
may  be  quite  difficult  to  construct;  hence  we  shall  assume 
that  I  p  is  not  known  and  that  the  first  term  on  the  right 
of  3° lo  must  be  bounded.   To  obtain  the  desired  bound  we 
introduce  the  auxiliary  function  T7   defined  by  2.39*   Since 
the  differential  operator  of  this  section  is  identical  to 
that  of  Theorem  2.2,  we  conclude  that   OC   and  Xr  may  be 
fixed  so  that 


3.50  <&  >   0 
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(x,y)£V 


3.51  J(9)<0  (x,y)£V 

and 

3-52  J(<VO<0     ^  >  O 

for  every  £>0  and  y  satisfying  0  <  £  *  y  -y  £  y  .   Next  we 
define  the  adjoint  J  of  J  by  writing 

3.53  J(  ?>)  =  J(^  )  -  Kg 

Because  the  singularity  of    7*   at  p  is  of  a  slightly 
lower  order  than  the  singularity  of  the  parametrix  y  ,  it 
follows  by  the  proof  of  Theorem  2.1  that 

=  0 

since  Z  is  strictly  positive  and  oC  is  a  non-negative  number. 
By  means  of  Green's  Theorem  we  write 

3.55  ju  m{^W)-^u^)}iY 

Wo) 


5k 

where  we  have  used  2>°5h   and  the  fact  that 


rr.i 


1<* 


^  ~z  *W 


-  #{*■&-?*  g} 4 


s 

for  p£D(y  ).   If  we  assume  that  K  is  non-negative,  then 
expansion  of  3»55  yields 

2  - 


3.5?  IU     JJJttfWCP)  -  V  7(9) 


+iMtm?h  -SffgJs  -  shifts 

S  S  -5  3 


+Sfr&z**f(i*<it)Js  =^ 


where  in  the  right  inequality  we  have  for  definiteness  set 
a  =  1  and 


S        V 

By  means   of   3*7    and   3°  21  we   have 


3.58  M6   =  max  {o9l+2M^F  -  ^  -rn^ ' %} 


3.59  itf2V?W)«lv  =   ZSSS^PF^plKU^)]^ 


Also 


3.60  £fJT<^l?rH-f>G<J  M,  W.7W) 


where  M.   is  given  by  3°  3^   and  OC  ^  is   an  arbitrary  positive 
3  P 

number.      Let 

3.61  cC  -  -2L 

5        2a 
o 

then  by  choosing  K  so   that 

M 


3.62  K^ja^  +  2iy=K2 

and  using  3»&0,  it  is  evident  that  3*57  may  be  rewritten  as 

3.63  fij*.  Iff  {*F  * fef  fcCfl.-^ -%)|^    <  Jt  • 


r&  v<4> 
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With  "9^  defined  on  D(y0)  -  {p}  as  is  indicated  in  Chapter  II 
and  with  the  constants  o ~  and  P  0   chosen  so  that  ft    *  P    and 
5  2^S0s  we  have  from  2.35  for  appropriately  defined  constants 
Co,  i  =  10,  .  .  .  13, 

3'6^      Jk*»     Iff  {&?Fj*fh(*ri)]-f?&)}<L)/ 
-  r(i*0]>^[^(v'V>]  -  P*  XflF)}  A 


=  0 


In  arriving  at  the  last  result  we  have  integrated  by  parts 
first  with  respect  to  ^   and  then  with  respect  to  y  (see  the 
proof  of  part  a  of  Theorem  2.1).   Prom  the  above  it  follows 
that 

3.65       /^?^[Vv^j(fJjll/i  y    . 


Prom  Schwarz's  inequality  and  the  fact  that  i1    satisfies  3-51 
it  follows  that 
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3-66  - \SfSW F  a^f,  [K(.%-*)]Jv\ 

>-{~SSS  vxJ(r)Jv}*- 


The  factors  on  the  right  of  3*66  each  exist.   The  first 
exists  since  t  is  continuous  on  V  and  is  hence  bounded  on  V 
while  1 (ty)    is  easily  seen  to  be  integrable  over  V.   For  the 
second  factor  we  note  that  the  singularity  of  *f   at  p€  D(yQ) 
is  less  than  that  of  J(  ^p)   and  the  result  may  thus  be 
obtained  direotly  from  Theorem  2.1.  Due  to  the  above 

3-67  {SSJ-3  (<d  ?*Jv}*<  { jJ/[- Wl',(?ffy$iyrj))lv]i 

and  then,    after  substitution  of  3«6£  and  3*20,   we  write 
3-68  {-J|{  7(<T)  ^Jv}k 
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Because  of  the  form  of  -fr    and  V  ,  we  use  Schwarz»s 
inequality  to  obtain 


3.69  \lim.  iSS  t>T(.xfi)J\/ 


{SiS  [-«*">]  T0t)*Jv}' 


where  the  first  factor  on  the  left  exists  due  to  Theorem  2.I4. 
and  the  second  factor  by  the  remark  following  inequality 
3.66.   It  is  clear  that  3*69  is  the  desired  bound  for  the 
first  term  on  the  right  of  3»l6. 

To  obtain  the  desired  bound  it  is  now  only  necessary 
to  bound  the  second  term  on  the  right  of  3«l6o   We  write 

3-70        \H^    l'l'j'yp{f(M)wJvw)-f(x,M  ifv<e)}M 


59 

In  the  last  inequality  the  first  factor  on  the  right  exists 
since  for  c  •  appropriately  defined 

3.7i       M«?<F''*>r%.*k-*)]J«  i    C,¥    t 

Is 

for  £  sufficiently  small.  The  integral  in  the  above 
expression  is  seen  to  exist  by  means  of  integration  by 
parts  as  in  the  proof  of  part  (b)  of  Theorem  2.1.   The 
last  factor  on  the  right  of  3. 70  is  now  bounded  by  recon- 
sidering 3-57  and  3»59»  Prom  these  inequalities  we  obtain 


3.72 


In  the  last  equation  the  choice  of  the  oC4.  'a,   depend  on  the 
magnitude  of  M^a   and  the  boundary  data.  For  definiteness 
we  set 

3.73  0^  *  1 

and 

3.73  a       <*2  .  5M3(ao-M3)"1      for  \>  VLf£ 


or 


3-73  b  oC     =  Ma"1  for  M  a_1>  i. 

d         j>  o  3  o       2 


Next  we   choose 


3'7l+  a  K  =  ab^-^L^^o-^1"1*2^2}38 


3 
or 


3.74  b  K  =  ~Jl+2i+Maal+a  M^M  2]=  K 

2nu  I         13o       o31J         • 


whence   either 


3.75  a  J^j  UfkR,z?9'x*zn3a;"£<iitp  ZAjv 


or 


3.75  b 


fes  #^>^iv;y^^^^i/ 
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where  S^   is  given  by  3-57  with  K  =  K~. 

Subject  to  the  remark  following  3*71  it  is  clear 
that  equation  3»70  and  the  proper  form  of  3*75  provide  a 
bound  for  the  second  term  on  the  right  of  3«l6. 

By  combining  the  results  of  this  section,  the  solution 
of  3°1  is  bounded  at  the  point  p£D(y  )  by 


61 


3.75-5       lwo„p-X(p)| 


*,./>* 


+1%  *  tijtw)}  'W*r  ['w-<t>]Jv}' 


In  the  last  inequality  the  functions  o   and  */»  are  given  by 
2.10  and  2.39,  and  the  constants  K±   by  3.I+2,  3.62,  and  3.7l±. 
The  computable  bounds  ~V>  ,  Q,   and  §  are  given  respectively 
by  3»1|8,  3-57,  and  3.75  while  the  function  X  (p)  is  defined 
by  3-17. 

If  problem  3-1  admits  a  solution  w  which  is  twice 
piecewise  continuously  diff erentiable  in  y  for  (x,y)6V, 
then  theoretically  the  pointwise  bound  may  be  made 
arbitrarily  small  by  improvement  of  the  approximation 
function  IP  .   In  this  regard  we  note  that  the  right  side 
of  3.75*5  is  zero  if  P  =  L  =  G  =  0.   Furthermore,  due  to 
the  form  of  the  bound,  it  is  clear  that  the  result  may  be 
Improved  by  means  of  the  Rayleigh-Ritz  technique  provided 
the  problem  is  linear. 
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By  referring  to  the  proofs  of  the  appropriate  theorems 
of  Chapter  II  and  the  methods  of  this  section,  one  observes 
that  the  boundary  value  problem  given  in  3*1  may  be 
generalized  by  allowing  the  components  a1J  =  a^J(x,y)  to 
be  non- symmetric.   Such  a  problem  was  not  considered  since, 
even  though  the  techniques  are  exactly  the  same,  the 
individual  inequality  expressions  are  naturally  much  more 
involved.   In  the  next  section,  where  the  matrix  a* J  is 
required  to  be  symmetric,  a  method  is  given  for  symmetrizing 
a  general  problem. 

As  a  second  generalization  it  is  easy  to  see  that  the 
differentiability  of  the  a1^,  i, j  m   1,  2  .  .  .  N  may  be 
reduced  by  merely  requiring  that  *»<   and  a*^  be  continuous 
in  the  interior  of  a  finite  number  of  regions  the  sum  of 
which  is  D(y)  for  o  <y<yQ»   The  a1  J  are  themselves  con- 
tinuous in  D(y)  for  o^y^y  and  piecewise  continuous  in  V. 
Under  this  reduced  differentiability,  one  may  easily  show 
that  the  techniques  of  this  section  are  valid. 

2.   The  Non-Linear  Dirichlet  Problem 

Our  aim  is  to  bound  the  solution  of  a  normally 
parabolic  problem  which  is  such  that  the  solution  is 
initially  given  and  thereafter  satisfies  a  non-linear 
normally  parabolic  differential  equation  in  the  interior 
while  assuming  given  values  of  the  boundary. 

We  assume  that  V  is  an  N  +  1  dimensional  domain  with 
boundary  D(y0)  U  S  t/ D(o)  which  is  such  that  2.1j.  and  2.5  are 
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satisfied.   The  boundary  value  problem  may  be  stated 
explicitly  as  follows: 

/J(w)  =  (a1^  w,1),1  -  Z-^2hf(x,y,w,vw)   (x,y)£V 
3.76       <w(x,o)  -  g(x)  (x,y)£D(o) 

^w(x,y)  =  h(x,y)  (x,y)€  5 

where  w  is  continuous  in  VUD(y0),  continuously  diff erentiable 
in  V U  S  and  is  such  that  its  second  derivative  in  x  and 
mixed  derivative  in  x  and  y  are  continuous  in  the  interior 
of  a  finite  number  of  subregions,  the  sum  of  which  is  D(y) 
for  every  y  satisfying  o<y<yQ.   In  addition  it  is  necessary 
to  assume  that 

film   w(x,y)  =  g(x) 
(x»y)-*(o,j) 

(x,y)  £V 


3.77 


lim  w,i(x,y)  -  g. (x)  i  =  1,  2  .  .  .  N. 

(x,y)— ^(ojy) 

U,y)£V 


The  differentiability  conditions  which  the  solution  of  the 
present  problem  is  to  satisfy  are  obviously  much  stronger 
than  those  which  were  required  for  the  solution  of  the 
boundary  value  problem  considered  in  the  first  section  of 
this  chapter.   Additional  conditions  must  also  be  placed  in 
the  coefficients  of  the  differential  operator;  in  particular, 
the  components  of  the  symmetric  matrix  ar*   =  a1J(x,y)  are 
piecewise  continuously  diff erentiable  in  V  and  satisfy  2.7 
while  the  function  2  -  Z(x,y)  is  piecewise  continuously 
diff  erentiable  in  y  for  (x,y)€V  and  satisfies 


3.78         min  2(x,y)  =  m?>o. 

V  c 

The  boundary  data  f ,  g,  and  g, .  are  assumed  to  be  integrable 
and  square  integrable  over  their  respective  domains  of 
definition  for  any  bounded  continuous  function  w  which  is 
continuously  diff erentiable  in  x  for  (x,y)6V.   Moreover, 
the  function  f  satisfies  a  Lipschitz  condition  in  all  but 
its  first  two  arguments •   Hence  there  exist  positive  numbers 
M.  and  M  ,  i  ■  1,  2  .  .  .  N,  such  that 

3»79        |f  (x,y,w1,  vw1)-f^x,y,w2,  v  w2)  |  <  M,  |w1-w2| +M1  |w1  ±-w2   A 

for    (x,y)  6  V,-  c^  <  w,vx  <  00. 

We  remark  that  the  problem  defined  above  and  the  additional 
problem  which  we  shall  obtain  from  it  by  introducing  an 
arbitrary  function  are  very  similar  to  a  problem  for  which 
Friedman,  [10]  and  [12] ,  and  Zeragiya,  [27]  and  [28] 
established  the  existence  of  a  solution. 

The  desired  bound  for  the  solution  function  will  be 
obtained  at  the  point  p6D(yQ).   As  in  the  previous  section 
we  assume  that  2  satisfies  2.9  at  p  (notice  that  2.8  is 
satisfied  due  to  the  differentiability  of  the  a1^)  and  that 
f  is  bounded  in  some  neighborhood  of  p  which  is  contained 
in  V  for  any  continuous  function  w  which  is  continuously 
differentiable  in  x  for  (x,y)£V,  that  is,  we  suppose  that 
3.5  is  satisfied.   To  obtain  the  desired  result  we  choose 
an  arbitrary  function  10  -I0(x,y)  which  is  3uch  that  J(  if  ) 
and  IP  (x,y)  approximate  f(x,y,  tf  ,  7  ^) ,  g(x)  and  h(x,y) 
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respectively.   In  order  that  the  divergence  theorem  may  be 
applied  In  what  follows  we  further  require  that  Lf>  be  twice 
piecewise  continuously  dlfferentiable  in  x  and  piecewise 
continuously  dlfferentiable  In  y  for  (x,y)6V.   Also  we 
require  that  the  LP,   be  piecewise  continuously  dlfferentiable 
in  y  for  (x,y) 6  V. 

Following  the  techniques  of  the  first  section  of  this 
chapter  we  define 

3-80         4>  (x,y)  =  w(x,y)-  ^(x,y) 

and  then  compute 

/j(^  )=f(x,y,w,  vw)-f(x,y,  tp  ,  7Y>)+F(x,y)    (x,y)  6  V 

3.81  <t  (x,o)=g(x)-  if  (x,o)  =  G(x)  (x,y)6D(o) 

l^(x,y)=h(x,y)-  q>(x,y)  =  H(x,y)  (x,y)  S  S 

where 

3. 82  FU,y)    =  f(x,y,  LP,  Vi^)-J(iP). 
For  later  use  we  note   that 

3.83  |jm|$M,  M     +M1  |^.|  +  |F(x,y)|  . 

With  lit)    and  X      defined  respectively  by  2.6  and  2.10 
we  note  that  3»l5  Is  valid  for  the  present  problem.   Sub- 
stitution of  3.8l,  3.82,  and  3*83  Into  that  equation  then 
yields  after  transposition 
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3^  lw<v*;)-zcfOl*|iU*.    tf-f  ^?(OAl 


UKWJf^'J 


where 


^         Me)  =  vk,o  -  r;  h  WJs  -(U.  Sfli  Mv 

+Ify0£G*/.s  -  JT/*  yp?w  <fs. 


The  volume  integral  on  the  right  of  3»85  is  proved  to 
exist  by  the  argument  associated  with  equations  3.18  and 
3ol9<»   Hence  w(x0,y0)  is  the  only  unknown  term  of  the 
left  side  of  3.8I|.  and  our  desired  error  bound  is  then 
clearly  available  as  soon  as  the  expressions  on  the  right 
of  3° ®h   are  bounded.   Since  the  unknown  volume  integrals 
will  be  bounded  in  terms  of  known  quantities  and  the 
integral  of  the  square  of  the  conormal  derivative  of  ^ 
over  S,  the  last  term  on  the  right  of  3 •  8^4-  is  considered 
first .   For  this  purpose  use  is  made  of  3«8l  to  rewrite 
2.113  as 
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3-86      ismfjis 


3      y  *     <r 


+  %  iff  f   tt  {i  |J  +  «j 3,  Wj  7w)-%, ^ 7(C)  t  fifct)}^ 


where  /{  is  given  by  2.11^.   Since  f  satisfies  3.79  we  may 
write 

3-87        J7/f  ^{«5^  +f(K,^w,7w)-f^f17.»)+Fl^}<|v 


f 

V 


.*.,  X* 


V       *  y 


where 


3.88        M  =  max  |z2  Slcr1)2^ 
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and  C(L_  is  any  positive  number.   In  the  last  inequality  the 
constants  M^_  and  IVU  are  given  by  3«25>  and  3-34  respectively. 
Hence  3«86  becomes 


3.89       fUttfAs  <^^3(c,^c,^^)fJU\^h 


Due  to  the  differentiability  of  the  solution  of  3«8ls  we 
may  write  by  means  of  the  divergence  theorem 

3-90      ffj  s*  ;r(^  i-  iff  i^»'  c*  *  A 


• 


Application  of  3*83  to  the  left  side  of  the  last  inequality 
yields 


3-91 


ii/ii5^^i^;/iia{fl^itflXiHpi}^ 


« I  iXTdit^  +  ft  Iff  *%  >  £i  $(%?Jv 


where  the  t   0  are  arbitrary  positive  numbers.   Subsequently, 
when  we  use  3»9^  ^°  rewrite  3"  90$  "the  arbitrary  constants 
will  for  definiteness  be  defined  as  follows: 
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3.92         £i=  £2=  £3=^ 

The  first  term  on  the  right  of  3. 90  is  bounded  by  noting  that 


where 


3*94         Mq  =  maxj  X- 


3a1  J  ^a^lZ 


v  *  ay  ^y 

Using  2.102  the  last  term  on  the  right  of  3.93  may  be 
rewritten  and  bounded  by 

3.95  -kSfo*\Us*-tiS~\($%-lrSU\(UfJU 

5  a  *  s         y  s       d 


s  ^       3 


where 


3.96         M  =  max  j-n"1^}, 


To  estimate  the  second  term  on  the  right  of  3. 91  we  first 
write  for  an  arbitrary  positive  number  c£p 
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3. 97  ff   4  ^5  *  T^(rv)</s  +  ^  iT(lg)Vs  I 

and  then y  almost  everywhere  on  Ss  define  the  vector 

3.98  T  -    (T1,  *    .  .  TNS  T7) 

where  T  is  tangent  to  the  surface  S,  is  orthogonal  to  the 
vector  (t  jo)s  given  by  2.97>  and  has  unit  Euclidean  length, 

I06  « 

3*99        X_  (T1)2  +  T^  =1 

We  denote  a  directional  derivative  in  the  T  direction  by 


tf-^TS^t 


Since  the  space  is  Euclidean, 


3.101        n  =  n^  for  every  i 


whence 


3.102  WTVV  +  *$"$+ ('-*{)*  St- 

Transposing  and  squaring  3*102  yields 
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hi 


o-^fCA)  •{<£*&*(-&#? 


-  ^  (p;)  M-  op*  ^  (t^rmj 


l  (<-<)( 


Since  2.5  is  satisfied  on  S, 


3.10k        1  -  n2>0 

7 


and  hence  from  3-103 


^  {S^fjs  *  m/0  f  (H)Vs  *m-$'($U 


♦tf*^M»--pjW<fc 


where 


3.106  M1Q   =  raax{2n^[aon(l-n^)]  "  }. 

S 

By  applying  the  above  results  to  J>. 90  we  have 
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v  V 

3M 


♦(%  ♦  ft)  J^  *  V** 


In  3.107,    let 


3.108  0Co  =  Vm 


2         v     10 

which  minimizes  the  coefficient  of  the  integral  of  the 
square  of  the  eonormal  derivative.   For  what  follows  we 


assume  that  (Mq+  2.VM  )  is  a  strictly  positive  number; 
when  this  is  not  the  case,  obvious  modifications  are  made 
to  obtain  the  same  result.   We  set 

3.109  04  =  ^2   (M  +  2  \n~  )  -1 


and  then  rewrite  3»89  as 
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J  a  v  ^   v 


If  the  last  terra  on  the  right  of  3*110  is  neglected,  then 
the  only  unknown  terms  on  the  right  of  this  inequality  are 
the  second  and  third.   To  estimate  these  unknown  terms  we 
make  the  transformation  of  variables  given  by  3*20  and 
observe  that  3»29  is  valid  for  our  present  problem.   In 
place  of  3*35  we  write 

J  a 
V 


3.H1  0-%g)  JT/.'^.M" 

0  \J  A      ( 


+(«mriMlt-Ml-&-%)sss!rUv 


<  «* 


*s 


jf/($)Vi  ihSWHffofardys-iSWtJs 


Ik 


+  H(  G^  °^P  (****>)  J s  +  *<x3  ^'-'^fg^-^v 


*  a 


In  the  last  inequality,  let 


3-.112 

3.113 
3- 111* 


3.H5 


a 
^2   M 


S"1 


J  M,   M,   -*«s 


<J#-T,{^  +  ^^HVJ^)*^+* 


r'/UJ3    j-HeST1 


^^^rm^m 


and  then  substitute  the  resulting  inequality  in  3«110« 
We  note  that  the  following  terms  appear  on  the  right  side 
of  3.110  after  this  substitution 
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3.n6  -^fe  +  ^V™*)*?** 

-  iU  (/%  +2  V7£)"'  //  * ' *  K;  ^  <ls 

Wo  '     m  *      * 


where  M,,  and  M  _  are  defined  by  the  above  inequality  and 
3.117        nL  =  min  Z(x,y)  >  o. 

Since  we  are  seeking  an  upper  bound  for  the  integral  of  the 
square  of  the  conormal  derivative  over  S,  it  is  clear  the 
unknown  terms  which  appear  in  3.II6  may  be  neglected  in  the 
revised  version  of  3. 110.   In  the  original  statement  of  the 
boundary  value  problem,  the  solution  was  assumed  to  be  con- 
tinuously differentiable  in  VUS;  if  we  make  the  stronger 
assumption  that  the  continuous  differentiability  holds  in 
VUSUD(y  ),  then  it  is  easy  to  obtain  an  upper  bound  for 
3.1l6  which  is  less  than  zero.   In  what  follows  the  stronger 
assumption  concerning  differentiability  is  made  since,,  if 
this  condition  does  not  apply,  then  one  may  clearly 
neglect  3*ll6  in  order  to  obtain  the  desired  result.   We 
define 
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3.118       B(yQ)  -  D(yQ)~D(yo) 


and  then  write 


where 


3.120        MD  =  max 

2  d(70) 


tfj. 


Next  we  set 

(M2+  C3Ch)  a   oC  -  M12 

3.121  — g      °      2      =     _ 


and  solve  for  c(  -  to  obtain 

5 


«,=±K'^&F). 


3.122 

Prom  the  above  it  is  clear  that 


3.123       -H,.  L~    N*    %,%!  Js  -A/.,  If  t  Js 
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The  inequality  on  the  right  of  3-123  is  strict  since  in 
Chapter  II  the  auxiliary  functions  f  were  chosen  so  that 
f*n.  had  a  positive  minimum  on  S. 

By  combining  the  above,  we  are  able  to  rewrite  3*110 
as  follows 

5  J   V 

-  i  H™\  {Tt)M  +  ^  If  Y  fa.*1  0-*ip}  '(H)^ 


8(30) 


^3 
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where  OC  »  is  given  by  3„ll5  and  the  computable  constant  \0 
is  such  that 

3-125        13-0 

when 

3.126        P  =  H  =  G  =  0. 

Inequality  3«12l|.  is  a  bound  for  the  unknown  factor 
of  the  last  term  on  the  right  of  3 •  81|. .   Bounds  for  the  two 
remaining  terms  on  the  right  of  that  inequality  may  now  be 
obtained  using  the  procedures  of  the  first  section  of  this 
chapter.   With  70  known  it  is  clear  that  3»57  becomes 

3.127  Lk  US  {&  ? U *)  -  * x  j  (v)  I 

j?cl»     "(1) 


1  z 
< 


ij^  H  ^[Hiu^ys  +13 


s 


+ 


The  right  side  of  3*12?  is  computable  and  hence  the  desired 
bounds  may  be  obtained  in  exactly  the  same  manner  as  was 
employed  in  section  1  of  this  chapter.   The  results  of  these 
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where  oC  .  is  given  by  3°ll5  ^.nd  the  computable  constant  (J 
is  such  that 

3.125       13  =  0 


when 


3.126        P  =  H  =  G  =  0. 

Inequality  3«12l|.  is  a  bound  for  the  unknown  factor 
of  the  last  term  on  the  right  of  3 .  8I4. .   Bounds  for  the  two 
remaining  terms  on  the  right  of  that  inequality  may  now  be 
obtained  using  the  procedures  of  the  first  section  of  this 
chapter.   With  10   known  it  is  clear  that  3*57  becomes 

3.127        Lk  ///  {&  *  n  *)  -  *  *  j  <*p) 


1    X 


JJT  W   ^f>[HlUV*)ys   *# 


^F ,  ,1 


I 


The  right  side  of  3*127  is  computable  and  hence  the  desired 
bounds  may  be  obtained  in  exactly  the  same  manner  as  was 
employed  in  section  1  of  this  chapter.   The  results  of  these 
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operations  may  be  written  as  follows: 
3.128        Iwfo^)-  X(p)| 


\ 


\  V 


'^WviS\^f^% 


IJTr-Tfrr^OVC^^*-^^ 


where  the  functions  Y      and  ^  are  given  by  2.10  and  2.39? 
the  constants  K-,,  IO^j  and  K  by  3.1lij.s  3„62  and  3- ?4 s. 
equation  3*127  defines  y^   and  the  constant  $  is  given  by 
3»75»   The  computable  bound  To   is  obtained  from  3 .124  while 
X (p)  is  defined  by  3»85«   As  in  the  problem  of  the  previous 
section,  we  remark  that  if  the  function  f  is  such  that 

3.129        f(x,y,w,vw)  =  f  (x,y) 

then  the  result  lends  itself  to  improvement  by  means  of  the 
Rayleigh-Ritz  technique. 

In  obtaining  the  conclusions  of  this  section,  the 
role  of  inequality  2.113  was  basic.  By  referring  to  the 
derivation  of  that  inequality  it  is  clear  that  the  symmetry 
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of  the  components  of  the  matrix  aiJ  was  essential;  hence  the 
consideration  of  a  method  for  the  symmetrization  of  an 
arbitrary  problem  is  appropriate.   The  result  is  obtained  by 
defining 

3.130  a1-*  ■  ^(a1^  +  a^1) 
and  then,  since 

3.131  aiJ^'ij  =  alJ^'ji 

almost  everywhere  for  (x,y)6V,  we  may  write 

3.132 


t 


where  3«132  is  also  valid  almost  everywhere  in  V.   Due  to 
the  foregoing,  the  differential  equation  of  the  boundary 
value  problem  may  be  replaced  almost  everywhere  by 

3,133        (5'JOr*^j r  ^^)  f%V*^ 
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Prom  3*133  It  Is  clear  that  ?  satisfies  the  same  conditions 
as  f  and  thus  the  original  problem  has  been  symmetrized  as 
required. 

Before  concluding  our  consideration  of  normally 
parabolic  problems,  we  remark  that  when  the  differential 
operator  has  the  form 

3.134  J(^)  =  (*1J  ffJ)fl  -^-ZyjF  ' 
the  transformation  given  by  3*20  yields 

3.135  J(?)  =  (arj  ?  )  t   -  (<fl  +  K)V     -Z^f* 

Hence,  if  the  minimum  of  the  function  as   over  V  is  a  strictly 
positive  number,  then  the  constants  K*  of  either  this  section 
or  section  one  are  reduced  by  the  amount  of  this  minimum. 
In  order  to  appreciate  the  Improvement  which  results  under 
such  conditions,  the  reader  is  referred  to  the  example 
presented  in  Chapter  VI. 


CHAPTER  IV 

DEGENERATE  PROBLEMS 

1.   The  Dir-ichlet  Problem  at  a  Parabolic  Point 

The  problems  which  have  been  considered  thus  far  have 
been  normally  parabolic  in  character.   We  now  turn  our 
attention  to  problems  which  are  normally  parabolic  at  certain 
points  and  degenerate  at  others.   As  the  above  heading 
implies,  this  section  pertains  to  the  derivation  of  a  bound 
for  the  solution  of  a  degenerate  boundary  value  problem  at 
a  normally  parabolic  point.   We  seek  to  bound  the  solution 
of 

/J(w)  =  (a1Jw?i)r  rqw-Zy|=f(x,y)       (x,y)6V 
1|-1        <W(x,o)=g(x)  (x,y)6D(o) 

\W(x,y)  =  h(x,y)  (x,y) 6  S 

where  w  is  continuously  diff erentiable  for  (x,y)  €  VU  S  U  D(y  ) 

and  is  such  that  its  second  derivative  in  x  and  mixed 

derivative  in  x  and  y  are  continuous  in  the  interior  of  a 

finite  number  of  subregions,  the  sum  of  which  is  D(y)  for 

every  y  satisfying  o<y<y  .   In  addition  It  is  necessary 

to  assume  that 

lim  w(x,y)  =  g(x) 

(x,y)-^(x,o) 

(x,y)€V 

lim  wfl(x,y)  =  g±  (x)  i  =  1,  2  .  .  .  N 

(x,y)— ^(x,o) 
(x,y) £  V 
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The  components  of  the  symmetric  matrix  a*J  satisfy  the 
conditions  Imposed  in  the  second  section  of  the  previous 
chapter.   The  functions  q  and  2  are  piecewise  continuously 
differentiable  in  y  for  (x,y)£V  and  Z  satisfies 

4-3         ZU,y)£o  (x,y)6V. 

Prior  to  stating  the  additional  conditions  which  are  to  be 
satisfied  by  the  functions  2  and  q,  it  is  necessary  to 
define  a  subregion  E  of  V  by 

k-k  E  ■  {(x,y)|  Z(x,y)  =0,  (x,y)6v}_ 

We  suppose  that  q  has  the  form 

where  the  constants  £  ,  I  =  1,  2  .  .  .  cc ,  satisfy 

&         r 

M       I  ll  -  I 

I'- I 

and  the  functions  q  ,  I  =  1,  2  .  .  .  it ,  fulfill  certain 
conditions  which  will  be  prescribed  in  what  follows.   Let 
the  coordinates  x*,  I  =  1,  2  .  .  .  N,  be  chosen  to  satisfy 

1|.7  min  max   ^(x1-**)   =  R 

(x»y)(x,yHS  l 

where   o  S  yl  y   .      For  N  =  1  we   suppose   that 

I|.8   a  m  =  min|o;q>  -aQ(l+R2)      } 

E  r 


8fc 


with  E  C7   defined  by 

o 


k>9       E  =j(xsy)   ||(x,y)-p||<  g   for  some  p6E  with  ^ojHV 

In  4-9  the  symbol  ||  ||  denotes  the  Euclidean  distance  in 
the  2  (and  later  N  +  1)  dimensional  space  of  the  problem. 
If  on  the  other  hand  N  =  2,  then  q  must  satisfy 

1^.8  b        m  ^injosr-j.q-^-aoClj.Rj}"1]   I  =  1,  2  .  ,  .& 
1  E 


4f(xi-4)2il 


where 

4.10  r 
and 

4.11  RT  =  max  r   , 

1    V   1 

Finally  for  N  2  3 

4.8  c        mjsminjojr^q  >-£  aQ(N-2)2  j    I  =  1,  2  .  .  .  X 

Since  the  techniques  which  are  used  to  obtain  solution 
bounds  for  N  =  1  or  N  =  2  are  very  similar  to  those  which 
are  employed  for  N-^3*  the  method  shall  be  presented  in 
complete  detail  only  for  the  latter  case.   Thus  for  N$3 
we  define  m  by 


4.12  os-ke1^ 


a0CN-2)2 


•1 

=  m  <1 


If  N  =  1  or  N  =  2,  similar  constants  are  introduced. 

Previously  we  required  that  Z  be  piecewise  continu- 
ously diff erentiable  in  y  for  (x,y)  6 V;  in  addition  to  this 


65 

it  is  also  necessary  that  -^§  be  continuous  in  V  for  every 
(x,y)  €  E  where  E  is  the  boundary  of  E„  and  §i  =  o  on    E/7(V-V). 

In  l|.l  the  boundary  funotion  f  «  f  (x,y)  is  continuous 
in  y  for  (x,y)€V  and  the  first  derivative  of  f  with  respect 
to  y  is  continuous  in  the  interior  of  a  finite  number  of 
subintervala,  the  sum  of  which  is  the  line  (x,y)  for  every 
fixed  x  where  (x,y)6  V.   In  addition  we  suppose  that  -r~ 
is  integrable  and  square  integrable  over  V  and  that  g  and 
g,,  are  integrable  and  square  integrable  over  D(o) .  By 
definition  it  is  clear  that  f  is  integrable  and  square 
integrable  over  D(o),  D(yQ),  S  and  V. 

It  is  now  our  aim  to  determine  a  bound  for  the 
solution  of  Ij..l.  As  the  heading  of  this  seotion  indicates, 
we  assume  that 

1^13        2(x,y)>o  (x,y)*D(y0) 

and  that  2.9  is  satisfied  at  the  point  p6  D(y  )  at  whioh  a 
bound  is  desired.  The  arbitrary  funotion  ^  =  if  (x,y)  is 
chosen  as  in  the  previous  chapter  and  the  function  f  (x,y) 
is  introduced  where 

Iplil        V(x,y)  =w(x,y)  -  tf(x,y) 

We  compute 

/j(Y  J  ■  P(x,y)  (x,y)£  V 

4.15       h   (x,o)  -  0(x)  (x,y)£  D(o) 

\?(x.y)  ■  H(x,y)  (x,y)f8 
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Let 

4-16         lit)    -  (a1^  )    =qt  +4(Z^' 

j        ^  y 

and  then  with  b   given  by  2.10  we  note  that  3»l5  is  valid 

P 

for  the  present  problem0   From  this  it  follows  that 

^17  Iwc^,^-^^)! 


i^  If!  ^  J(8f)  J^ 


^{IK^ilMl/yAf 


where 


fc.ie  x(f)«  ^0^0~^!f^ 


-fium.  ill  ^Fjv  +  If  Y*2G^  -  ttnXzHds  , 


Since  P  is  bounded  over  V,  it  is  clear  that  the  volume 
integral  In  L..1.8  exists 5  thus  the  desired  bound  is  available 
when  the  terras  on  the  right  of  I4..I7  are  estimated.   We  shall 
first  consider  the  integral  of  the  square  of  the  conormal 
derivative  over  S  and  then  use  this  bound  to  estimate  the 
unknown  volume  integral.   For  this  purpose  we  make  the 
transformation  of  variables  given  by  3»20u   The  boundary 
value  problem  then  becomes 


8? 


4.19        ?(*  0)  =  g  ^  ( n^  (y^eD(o) 


where 

4.20         q=KZH, 

^  Z  - 

Due  to  the  fact  that  -*-=  is  continuous  in  V  for  (x,y)€  E 

o  y 

and  I4..8  is  satisfied,  we  may  choose  K>o  so  that 

for  N  ^  3  and 

4.22         min  (±KZ  -  ±  4-5)  >  -  (1-5)  -° 

V    2      2  ^y  ^Z 

where 

14--23         5  =  KZ+f  XqT=  Z_£I(Jr2f+q-r)::=£1^ 

In  the  above  we  have  made  use  of  the  fact  that  -S-=  =  0 
on  E.   Similarly  the  constant  K>o  is  chosen  so  that 

l[.21  a       m  =  minfq-§ZK>-a  (l^R2)"  }o(   for  N  =  1 
and 
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4.21b  ^r  =  />r^{o     '  ^    (f    -A*^LrJ"')>-^(V^)"/] 

With  K>o   as    chosen  above   we  have   from  2.113 


V 
V 


+ ^'  jut*  ( ^^/  Jv 


where  dt  -  is  an  arbitrary  positive  number  and 
I4.25         M13  =  maxj"z  Ijf1)2] 
L,  .26         M,  1  -  max  |q  j 

In  order  to  estimate  the  unknown  terms  on  the  right  of  l±»2l\. 
we  observe  that  for  any  set  of  N  functions  g1,  i  =  1,  2  .  . 
.  Ns  which  are  piecewise  continuously  diff erentiable  in  x 
for  (x,y)6v,  it  follows  that 

By  setting 

1+.28         g1  =  (x^x1) 
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In  4.27 ,  it  is  easily  shown  that 

Substitution  of  this  inequality  into  1+.2I+  yields 

4.30         ucdfds 

To  bound  the  unknown  terms  on  the  right  of  ij..30  we  write 


-///?  ^W^)]^  tajK^^^-^Ww 


+f  sf  $*  -*il$  V°Js  -i^fHV^MV- 
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Did) 


Due  to  the  assumptions  which  were  made  concerning  the 
differentiability  of  the  solution  of  l±.l9    it  is  clear  that 
3° 93s  3°95>  3 -97s  and  3-105  are  valid  for  the  problem 
presently  under  investigation*   Substitution  of  these 
results  into  a  transposed  form  of  I4. «, 31  gives  after  making 
use  of  the  fact  that  ^  =  P    on  D(y  ) 


^VV4^^W'A 


where  we  assume  that  (Mq  +  2  VmTq)  is  a  strictly  positive 
number.   If  this  is  not  the  cases  obvious  modifications  in 
what  follows  will  yield  the  same  result.   In  I4. -  32  it  is 
clear  that 
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and 


+ 

y 


is 


1+  •  3L.         M,  *  -  max  V  4^ 


We  now  show  that  due  to  the  form  of  q  and  henoe  q  It  la 
possible  to  choose  £.  so  that  the  terms  in  the  braces  on 
the  left  of  4* 32  may  ka  bounded  from  below  by  a  computable 
constant.   On  the  surface  D(y0)  we  have  by  means  of  the 
divergence  theorem 
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-  If  %  Ms  +z(S  V  %  <?M 


where  B  (y  )  is  the  boundary  of  D(y  }.   For  N^3  the  form  of 
o  o 

1|.35>   leads  us    to   define 

4.36  g1   =    (x1-xi)r-| 

and  then  for  every  I  obtain  from  1^.35 


where  cC  is  an  arbitrary  positive  number  satisfying 

We  choose  cC  so  that  the  coefficient  of  the  last  term  on 
the  right  of  I4. . 37  is  minimized.  The  above  equation  thus 
becomes 
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4*38  o 


*&  JV-#«?  *,  //?  *  fa-  SI  a.'  ** «  M . 


For  later  use  we  observe  that  It  Is  also  possible  to  show  that 


k  -39  o  fffP^Uv 


fattttyf       ** 


^^+<foaW*X^ 


For  N  =  1  we  have 


nx**, 


and 


4-39  a  ff  ?Vv    <   J  /(«  -*K  ^   *  *£  #«"  «  %  JV 

V  i  *  He     V  *     * 

Finally  for  N  =  2  we  may  write 
and 

4.39  *      mtyitttfwy^Mi+frjfft'i.M* 


9k 

We  are  now  in  a  position  to  obtain  a  lower  bound  for  the 
terms  in  the  braces  on  the  left  of  4.32.   The  details  shall 
be  worked  out  only  for  N^3«   Similar  results  easily  follow 
in  the  other  two  cases.   By  means  of  the  above  and  \\.2\   c 


We  notice  that  the  line  integral  on  the  right  of  l+.I).0  is 
actually  a  volume  integral  over  an  (N-l)  dimensional 
domain;  and,  since  N£3,  it  is  clear  that  the  integral 
exists.   In  i|.32  we  choose 

N2a, 


4.1*1         E1   =  (l-m)-j-|>0 


where  the  strict  inequality  on  the  right  follows  from  l\..12, 
Substitution  of  I4. •  I4.X  into  \\»3>2   produces  the  inequality 

v   V*   -.   -y 


w*        m^h  *UH,«*%.)*p(**$ss$fi 
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«*) 


In  deriving  the  above,  an  inequality  analogous  to  I4.29  over 
D(y0)  is  used  to  bound  the  integral  of  the  square  of  ^  over 

D(y0). 

So  that  I4.3O  may  yield  a  non-trivial  bound  for  the 
integral  of  the  square  of  the  conormal  derivative,  we  choose 
0C-,  as  follows: 

U-U3  #,  •-|m3(M9+2VrF'o)''1exp  (-2KyQ) 

Inequality  14.30   then  gives 

3        y  *      '$  K- 


where 
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■j   i 


The  unknown  term  appearing  on  the  right  of  l±*l\k    is  estimated 
by  using  the  divergence  theorem  to  write 

V  r  V 


+fi 


i/;  f  VN»-# » +  ^  #  <■ Wfe 


~££r^  (&)  f  (**-*i)^SH  WN^y5 


+^  X0  a.    £.  ^  </v    f-fc  |T  £  fr  >^  (*  a  ^)  ^s 
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*k.r, 


where  L  -,  and  £   are  arbitrary  positive  numbers  which  shall 
be  suitably  chosen  in  v/hat  follows.   To  obtain  the  above 
it  is  clear  that  I4..21  c,  Ij..22s  i|.29  and  !|o38  c  have  been 
used  as  necessary.   Inequality  L[.l|6  is  now  transposed  so 
that  the  unknown  term  on  the  right  of  \\cl\l\   may  be  bounded 

k'kl  *(<-*)  iff*"?, .*Am 

r-i        1        '  s 


£  /^.  (A/VJj(^'xrJ^VHWp[^K^0-^)]^o* 


id 
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where  n)  i  Is  defined  by  the  above  equation  and  It  is  clear 
that  we  have  set 

1      a      R2 

Since  in  <  1  due  to  If. 12,  Inequality  if.lf7  provides  a  non- 
trivial  bound  for  the  Integral  of  a  ^  ^   ±   t  *•      Bv  comparing 
If.lflf  and  h'kt   we  are  lQd  to  choose 

then  /f.lfif  becomes 

where  7j  la  computable  and  is  such  that 

If. 5i       13  -  o 
if 

If. 52         P  «  H  *  0  »  0 

The  above  Is  a  bound  for  the  unknown  factor  of  the 
second  term  on  the  right  of  If •  17 •  W©  next  bound  the  volume 
integral  on  the  right  of  that  Inequality  In  terms  of 
computable  constants  and  If .50.  Let 

If. S3        J*(lP]  ■  (alj  ^#1)   -  2q^  +^  (a^) 
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and  consider  the  function  T  given  by  2.7l{.*   Since  the  lower 
order  terms  of  J^  and  "J  are  identical,  we  may  deduce  from 
Theorem  2.3  that  the  arbitrary  constants  of  2.7lj.  may  be  chosen 
so  that 

k.5k  -PiM      and   J*(<f)£o       for  (x,y)£V 

and 

for  every  (x,y)€  V  such, that  o<l<   y  -y  i  y  .   Prom  the 
proof  of  Theorem  2.1  and  the  definition  of  f    and  2V  it 
follows  that 

By  means  of  the  divergence  theorem  and  the  differentiability 
of  n7    we  have 

t-57  JU  W  fW)-*VCF)}  A  -10*  t^Jv 


Wo) 
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and  by  expansion  of  the  last  equation  the  following  inequality 
is  established 


where  dJ  is  obtained  from  I4. •  I4. Y"  after  application  of  1|..5>0 
and  M  is  given  hy  \±,%l\.o  By  an  argument  similar  to  that 
which  was  used  to  establish  3.6I4,  we  are  able  to  write 

and  finally  from  equations  analogous  to  3»o6s  3*67*  3«o8 
and  3»°9  we  have 


*  {  kt  tcr^r  j  (&)%M\{  xr/[-T^)K^J  ^ 


^{^^[-t^)]V^^^ 


The  existence  of  the  Integrals  on  the  right  of  [(..60  is 
established  by  a  discussion  similar  to  that  of  the  previous 
chapter. 
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Inequalities  1\.»%1  and  I4. . 60  give  computable  bounds 
for  the  unknown  terms  on  the  right  of  1j..17  and  hence  the 
desired  pointwise  bound  is  obtained  from  the  inequality 


^•6l        Mv^-^Kifjr-jW'y^v] 


mm®? 


where  the  functions  o   and  "h   are  given  by  2.10  and  2.7ij. 
and  the  computable  constants  10   and  ^  are  obtained  from 
4.50  and  l|.o5>8.   As  in  the  previous  problems,  we  observe 
that  the  right  side  of  If.6l  is  identically  equal  to  zero  if 
P  =  G  =  H  =  0. 

During  the  investigation  of  the  problem  of  this 
section,  an  attempt  was  made  to  generalize  the  method  by 

omitting  the  requirement  that  -^ —  be  continuous  in  V  for 

7  ^y 

(x,y) 6  E.   Such  a  generalization  has  not  as  yet  been  achieved; 

however,  if  the  domain  V  is  restricted  as  indicated  below, 

then  the  desired  bound  is  available. 

We  suppose  that  the  given  boundary  value  problem  is 

identical  to  that  stated  at  the  beginning  of  this  except 

that  -4-=  need  not  be  continuous  in  V  for  (x,y)eE.   In 
«  7 

addition,  however,  we  suppose  that  V  is  a  cylinder  (that 
is,  D(o)  =  D(y)  for  o£y£y  ).   It  is  evident  that  equations 
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I^.llj.  to  1|.20  are  valid.   We  choose  K>o  so  that 

Z4-.62         mjamln-Tcr^q  >  -raQ(N-2)2j   I  =  1,  2  .  .  .^ 

for  N£  3  where  it  is  clear  that  K  >  o  may  also  be  similarly 
defined  for  N  =  1  and  N  -  2.   In  what  follows  we  shall 
consider  in  detail  only  the  case  N£  3  since  similar 
techniques  would  give  the  same  results  in  the  other  two 
cases. 

The  integral  of  the  square  of  the  conormal  derivative 
of  y    is  obtained  by  bounding  the  unknown  expressions  on 
the  right  of  \\*J>0   in  terms  of  the  given  boundary  data.   Prom 
the  divergence  theorem  and  the  differentiability  of  ifJ   and 
the  coefficients  of  the  differential  operator,  we  have  for 
an  arbitrary  constant  "a" 

^3       §!%*(*)   *f^r*)]^ 
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*(-%  +  & 


«l)  $«■ '  i  \  *f  fra>-#'  v-iSf  *X  fy  Js 


^X^h^^s  *  ^srcijjs 


~7 


If/If  f-l»f  |^.-a)]A  -  ^ *(s|)*f  ^^^i/. 


+t 


The  above  inequality  is  transposed  and  use  is  made  of  if. 29: 
if.  39  c,  if.  62  and  SchwarzTs  inequality  to  obtain 


if.6if 


Ja     tfft     a  ;r     x.^irKx    'Latere   « j  ~    ~    1 

**-s{-** -■*£■-  niJfa  ^^^ 


+ 


5 


^^  ^^^P^+M^^  4  "4^1  ^[^K)&]<4 


J)(.) 


.  JL 


lOij. 


In  1^. . 30   and  the   last  Inequality,    we   set 

MS 


m- 


1        2 


4.66 

4-67 

and 

4.68 


£x-i 


N2a, 


^2  -  J^F  C1  -  -) 


Substitution  of  the  above  into  4.3O  yields 


4.69 


tf  (  &S    *  1%  f  ^  XI(,  V)^,  *ty  gK  <^)]  * 


;i  ^fVN^v +xr[^-n#]*Y>t*^^.^)yi 


'W 


fiiiri?- 


*3-(«hof]  ifr^^-^]^  =  ■# 


where  a  is  given  by  I4. . 68  and  K>o  is  chosen  so  that  1|.62 
is  satisfied.   Since  Tj  is  a  computable  constant,  the 
desired  bound  for  the  integral  of  the  square  of  the  conormal 
derivative  has  been  obtained.   A  bound  for  the  unknown 
volume  integral  in  1|.17  is  again  given  by  4.60,  and  thus  the 
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desired  pointwise  bound  is  available. 

2.   The  Dirichlet  Problem  at  a  Degenerate  Point 

We  seek  to  obtain  a  bound  for  the  solution  of  i|«l  at 
a  degenerate  point?  that  is,  at  a  point  p  which  is  such  that 

h-10  Z(p)  -  0 

At  a  degenerate  point  the  author  was  unable  to  determine  a 
parametrix  with  properties  analogous  to  those  of  ^   at  a 
normally  parabolic  point.   Because  of  this  reason,  additional 
conditions  had  to  be  placed  on  the  boundary  data,  the  shape 
of  the  domain  and  the  region  of  validity  of  the  differential 
operator. 

Let  the  problem  for  which  a  bound  for  the  solution  is 
desired  be  that  which  is  proposed  at  the  beginning  of  the 
last  section  (specifically,  from  the  heading  to  the  paragraph 
immediately  preceding  equation  ij..l3)  with  4*70  replacing  Ij..l3< 
In  addition  we  suppose  that  V  is  a  cylindrical  domain,   Jv 
is  twice  piecewise  continuously  diff erentiable  in  y,  and  the 
differential  equation  is  satisfied  on  D(o) .   Moreover  we 
assume  that  the  derivative  with  respect  to  y  of  the  conormal 
derivative  of  w  is  continuous  in  the  interior  of  a  finite 
number  of  regions  the  sum  of  which  is  S  and  that  the  fol- 
lowing are  satisfied 

k-n  &n       *($>)-- IU       1(a) 


f,  6D(o)  p%es 
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If. 72        Z(x,o)  >  o 

4.73        llm-i^<P)  -H^W) 

37       dy 

£->(x,o) 

p  £V 

and 

lf.7^.         2(x,yQ)  s  o 

for  some  open  set  D'CD(y  )  for  which  p  is  an  interior  point. 

To  shorten  the  argument ,  we  consider  only  the  case 
N^3  since  for  N  =  1  or  N  =  2  similar  results  are  easily 
obtained.   As  in  the  previous  sections ,  we  obtain  the 
desired  bound  by  introducing  a  function  tf   =  If   (x,y)  with 
the  same  differentiability  as  w  and  then  with  ^    given  by 
If.llf  we  compute  the  system  !f.l£.   Next  we  define  the  function 
f    given  by  3.20  and  compute  If. 19*   The  constant  K  ^o  is 
chosen  so  that  If.  21  c  and 

1       1^5*  a0N 

If. 75         min  (~KZ  +  ~  -r-5)  >  -(1-m)— 5 
V   2     2  ^         8R 

are  satisfied.   The  desired  bound  is  now  obtained  by  using 
(6.8)  of  [2l]  to  write 


+  f  {^f-nyi}ji 


BQd 


IF 
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where 

k-n  &(t)  =   (VV,.),..  -^f 

and  P  is  given  by  (I4..2)  of  |_21|  .  The  unknown  terms  on 
the  right  of  \\»1%  are  the  first,  second  and  fourth.  For 
the  fourth  term  consider  that 

Using  Hormander's  Identity  and  the  same  techniques  as  were 
used  to  obtain  2.113,  it  is  clear  that 

3  %)  * 


where 


4.80        t'bihfWfct'ti-gpte 


_        -lk 
lj.8l        m„  =  min  n  f  n, 

3   B(y0)     k 


and  f   =  f  (x)  are  any  set  of  N  piecewise  continuously 


10b 


differentiate  functions  which  are  such  that  f  n,  is 
bounded  and  has  a  positive  minimum  on  B(y  )»   Prom  If. 32 
after  application  of  lf.f>0  and  If. If  7  or  from  lf.6If  with  the 
unknown  term  on  the  right  bounded  by  If. 69,  it  is  evident 
that  one  may  obtain 

4.82  ft    a'V  %  Is  £  D 

where  Js  is  a  computable  constant.  Inequality  If. 82  bounds 
the  last  term  on  the  right  of  if.?9>  for  the  second  term  we 
wri  te 

+S$*(ifJs  +  XT  f*^ 


where 


r^   «  2 


i+  •  8I4.         1^  =  max  I^L^f1)  j 


»(y0r 


If. 85         M_   -max  (zZc^1)2} 


dc^o) 


and 


109 


h . 86         M.  ,  =  max  {  q2  \ 


By  using  an  inequality  over  D(y  }  which  is  analogous  to  1^.29; 
we  may  rewrite  l\..79   as 


rr   o  1  .   — 


3  o  o  n   a0  3      B(^ 


For  the  last  term  on  the  right  we  use  the  divergence  theorem 
as  follows: 


It.  88 


Since 
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it  follows  that 


Prom  hf\\l   after  appl  I  caption  of  4»50  or  from  I4.  •  6I4.  with  the 
unknown  term  bounded  by  4»°9>  one  has 

By  transposing  I4. •  88 ,  and  integrating  by  parts,  we  obtain 


+  £lf(x -*">], I"**  ♦  S$]  *J»|fcKfy-a)]  <^ 


Sl  £fJ"[-*F+  5^J  B^p t?K(fy-^]  </v     +  % 


where 


Ill 


I--/ 


ll\   (a7t)//(xa-^\\{-^+  ^}^[^f#/i 


In  I4..92  we  choose  the  strictly  positive  numbers  CC ..  and  Q£ 
so  that  the  coefficient  of  the  second  term  on  the  left  is 
zero.   By  comparing  I4.87  and  1|.92  we  are  led  to  define 

By  fixing  the  arbitrary  positive  constants  as  indicated 
above,  it  is  clear  that  one  easily  obtains 

k-91  f  {¥vfJH   <   ft 

where  p  is  a  computable  constant.   Inequalities  1^.95  and 
Ij. . 78  provide  the  desired  estimate  for  the  last  term  on  the 
right  of  4 .75-   The  second  term  on  the  right  of  that 
equation  is  bounded  by  using  exactly  the  same  procedures  as 
Payne  and  Weinberger  employed  in  [2lJ  .   Let  r(p)  be  defined 
by  (I1-.I8)  of  [2l]  ;  it  is  then  clear  that   /O 


and  Xr 

0 
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may  be  chosen  so  that 


4.96 


and 


4-97 


-Zf  i  M 


for  (x,y) 6  D(yQ) 


(L*FF)     <0 


for  (xsy)£  D(yo) 


where 


4.98 


*(?)-(*%),<-*$¥ 


Prom  the  divergence  theorem  and  the  definition  of  ~r 


1+.99 


z~+ 


.2  3=F 


and  thus 


4.100 


^0) 


•k?A 
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To  bound  the  unknown  term  in  the  above  expression  we  make 
use  of  I|.82  as  follows 

ill1 

and  then  obtain 


4.102  tftftfrft'^actfjis 


where 


4  •  103  M_  A   =  max      {  t-a*  WC  )J  °VM 

4.104  M        =  max      {^ZK2} 

and  from  4»92  after  application  of  4*95 

4.105       jr  e  ( U)V^  3». 

By  means  of  Schwarz's  inequality 

4.106  \S$vatrp)Js\ 

In  the  above  the  first  factor  on  the  right  is  bounded  by 
4.100  and  4*102  while  the  second  may  be  computed  due  to 


ii4 

(/j.l  b)  of  [2l]  .   To  estimate  the  first  terra  on  the  right  of 
4.75  we  consider 

h-w  U  PM)  is 


Ditii 


where 


4.108  Jifl   =     max      (gp2) 

18     D(y0)-D*       PJ' 

Transposing  I|..75  and  taking  the  absolute  value  yields 


where  l|.  106  and  I4 •  78  bound  the  last  two  expressions  in  terms 
of  computable  constants. 

It  is  the  belief  of  the  author  that  the  conditions 
imposed  on  the  boundary  value  problem  of  this  section  are 
excessive.   In  future  investigations,  efforts  shall  be 
made  to  produce  pointwise  bounds  at  a  degenerate  point 
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without  Imposing  the  rather  restrictive  conditions  which 
were  necessary  in  order  that  the  above  techniques  yield  the 
desired  result. 


3.   The  Mixed  Problem  at  a  Parabolic  Point 

Thus  far  in  this  chapter  we  have  always  required  that 
the  matrix  a1^  be  positive  definite  and  that  S  be  a  non- 
negative  function.   Now  we  suppose  that  2  is  strictly  positive 
but  reduce  the  conditions  on  the  matrix  by  requiring  that  a  J 
be  positive  semi-definite. 

Let  V  be  an  (N  +  1)  dimensional  domain  and  let  the 
notation  of  the  problem  be  as  indicated  in  equations  2.1  to 
2.[|..   In  place  of  condition  2.5  we  now  require  that 


1^.110 


min  in—   -  m^>  o 


We  shall  derive  pointwise  bounds  for  the  solution  of 


4.111 


V 

where   clearly  it   is   necessary   that 


4.112 


H  (x,y)    =  0 


(x,y)  £V 

U,yK  D(o) 
(x,y)€  S 


for    (x,y) £  S±  with 
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If. 113        S1  -  {  (x,y)  |  (x,y)  €  S,  l_ni  =  °}  . 

The  solution  function  w,  the  coefficients  of  the  differential 
operator,  and  the  boundary  data  have  the  same  differenti- 
ability and  integrability  as  the  analogous  functions  of  the 
mixed  problem  of  Chapter  III.   As  in  the  previous  problem,  the 
function  f  satisfies  a  Lipschitz  condition  in  w,  that  is. 


t[.ll4  |f  (x,y,wx)-f  (x,y,w2J|<EL|w1   -  w2|     t 

As  the  title  implies,  the  solution  is  to  be  bounded  at  a 
parabolic  point  p6D(yQ}.   Thus  a1^  is  positive  definite  so 
that  a..(p)  exists.   In  addition  to  the  above  it  is 
necessary  to  assume  that  conditions  2.8,  2.9,  and  3*5  are 
satisfied  at  p. 

We  use  the  notation  and  techniques  of  the  previous 
mixed  problem  to  obtain  the  following  equation  which  is 
analogous  to  3»1° 

S  S  Q 


where 
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4.116  *(?)  "W <*>*,)*■  SfLV-  Js 

D(d)  j  y^1>     v(t) 

The  right  side  of  l±.ll6  is  computable  by  the  arguments 
associated  with  equations  3-l8  an^  3*195  an<3  hence  the 
desired  bound  is  available  when  the  unknown  expressions  on 
the  right  of  1|.  115>  have  been  estimated.   To  estimate  the 
boundary  integral  the  change  of  variables  given  by  3«20  to 
3.22  is  made  and  then  by  means  of  the  divergence  theorem  it 
is  clear  that 

h-  w  /r  if1  j  s  <  2  („*)"  {%SS9XM 


+ 


^.  ^  iZ  ^t  p-  K^-*)j J<* + (~h +«,+£+ i  "*)  Sit  Plfi/ 


^///^^[z/c^-^A+iX/^^Wf^^V^ 


where 


lj..ll8  min  2  -  nu^nu>o 

S  ^        ^ 

In  order  to  obtain  a  non-trivial  bound  for  the  left  side 
of  I|..117  we  set 


it.  119         \   =1(2^) 


and 
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U-120       /C«/m;'(/T,  +iHMv)M; 


Dae  to  the  above  we  may  write 


fc-121      K  <fjs    <4(<n,<msfSf?**f>[*W-ii-'$Js 


In  order  to  bound  the  volume  Integrals  of  I|.ll5  we 
introduce  the  function  -h   given  by  2.39«   Even  though  the 
differential  operator  of  this  section  does  not  satisfy  the 
conditions  of  Theorem  2.2,  an  inspection  of  the  proof 
reveals  that  the  conclusions  of  the  theorem  are  valid  for 
the  present  problem.   First  we  observe  that  constants 
f>    t      &£  ,  and  S<ft  associated  with  equations  2.36  and 
2.37  exist  since  J  is  normally  parabolic  at  the  point  p 
and  thus  a1 J  is  positive  definite  in  a  small  neighborhood 
of  p  due  to  2.8.   Because  2  is  strictly  positive  in  V 
it  easily  follows  by  remarks  similar  to  those  which  were 
employed  In  the  proof  of  Theorem  2.2  that  the  constants 
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or   and  xr   of  2.39  may  be    chosen   so   that 

1|.122  l&)io,    ^io  (x,y)£V 

and 

4.123  1(9)  <  o9^<0  ,  (x,y)e{x,y  |(x,y)ev,o<  ££  yQ-y] 

Due  to  the  foregoing,  equations  3.I4.9  to  3*75  ar-e 
valid  for  the  present  problem  with 

k .  124        M2  =  M1  =  t)  . 

The  bounds  for  the  unknown  volume  integrals  of  1j..115>  are 
hence  available  by  means  of  the  appropriate  equations  of 
Chapter  III.   After  application  of  these  results  one 
obtains 


ij.l2£        I *<*,*;  ~  ^Cf>)l 


v  r 
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where  a   p  is  given  by  2.10  and  "^  is  defined  by  2.39  subject 
to  the  remarks  of  this  section.   Using  the  results  of 
Chapter  III  we  obtain  that  Kp   and  K->  are  given  by  3,62  and 
3.7^a  while  the  computable  constants   y2  and  ^  are  defined 
by  3«5>7  a^  3«75a.   The  bound  "P  is  determined  by  the  methods 
of  this  section  and  is  given  explicitly  by  1|,121<> 


CHAPTER  V 
THE  ELLIPTIC  PROBLEM 

1.   The  Neumann  Problem 

In  [21^ 9   Payne  and  Weinberger  obtained  pointwise 
bounds  for  the  solution  of  second  order  interior  elliptic 
differential  equations  with  Dirichlet  or  mixed  boundary- 
data.   The  occurrence  of  the  zero  eigenvalue  of  the  free 
membrane  problem  prevented  the  authors  from  extending  their 
results  to  Neumann  boundary  value  problems.   By  following 
closely  the  procedures  which  were  used  to  obtain  a  bound 
for  the  solution  of  a  normally  parabolic  mixed  problems,, 
we  obtain  such  bounds  for  an  elliptic  problem  of  the  same 
general  type. 

In  what  follows  we  shall  use  the  notation  and  many 
of  the  results  of  section  6  of  [21J  .   We  seek  pointwise 
bounds  for  the  solution  of 

B(w)  =  (a1Jw  , )  .-b  w   -qw  =  f(x)        x£D 


S»  ■*■  S 


i 


5-1 


^W=  ^ 


where  D  is  a  domain  of  Euclidean  (V-space  with  boundary  B. 
In  order  that  the  divergence  theorem  be  applicable  it  is 
necessary  to  assume  that  w  is  continuous  in  D,  continuously 
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diff erentiable  in  D  and  Is  such  that  its  second  derivative 
is  continuous  in  the  interior  of  a  finite  number  of 
subregion,  the  sum  of  which  is  D.   In  addition  it  is 
necessary  to  assume  that  the  solution  of  the  boundary  value 
problem  is  such  that 

5.2  lim  aij°(p)w   (p)  =  a1  J  (x)w  , .  (x)    1b  1,  2  ...  B 

p_-»x  €B     J  J 

P  £D 

We  suppose  that  the  components  of  the  matrix  a  **  and  the  at 
are  piecewise  continuously  diff erentiable  in  D  while  q(x) 
Is  bounded  and  strictly  positive  in  the  same  domain.   As  in 
(6.3)  of  [2l]  it  is  necessary  that 


s-3  ^  k[Ti^(jig;U  ^ej^u'fj^^i 


a,[z^  +  1w^] 


for  all  real  numbers  ( £   ...  ?  )  where  equality  holds  on 
the  left  if  and  only  if  all  the  numbers  are  zero.   Let 

5.I4.         min  q(x)  =  m,  >  o 
D 


where  iil.  is  positive  since  q(x)  is  strictly  positive  by 
assumption.   From  (6.8)  of  [2l]  we  have 


D 


6 


123 

where 

5.6  B(  Pp)  =  (aiJ  f   )    +  (  J/,1  p   )      -q  PP 

Ps  J  » i         p  j  i 

with  Pp  given  by  (I].. 2)  of  [2l]  and  if   defined  as  in  the 

second  paragraph  of  section  l\.   of  that  reference,,   We  apply 

Schwarz's  inequality  to  the  first  boundary  term  in  £.5  and 

obtain 

where,  as  in  the  previous  parts  of  this  thesis,  we  have  set 

Let  f  (x)  .  .  .  f  (x)  be  the  set  of  functions  introduced  in 
sections  2  and  3  of  the  work  by  Payne  and  Weinberger,  that 
is,  the  set  f  (x)  is  piecewise  continuously  diff erentiable 
in  D  and  is  such  that  f  n,     is  bounded  and  has  a  positive 
minimum  on  B.   For  later  use,  we  set 

5.9         mP  =  min  fkn 

B 

and  then  bound  the  unknown  factor  on  the  right  of  5*7  by 
writing 

5.10         Ms  <-4 1  H^h  --i  j/Wl  J  v 

*  j)     *  p 


% 


12lf 


where 


5.11  M,    =  max  If1,  \ 

5.12  M2  =  maxJXf^J 


D 

and  £  1  is  an  arbitrary  positive  number.   By  means  of  the 
boundary  value  problem  and  the  divergence  theorem,  it  is 
evident  that 


5.13  JS<PfJv  =-N%fi**,Jl\,  +&UJ$ 


D  A  0     • 

where 

5.11+         F(x)  =  f(x)-B(</>) 

and 

5.15       L(x)  =^U)  -  -J7"  • 

By  means  of  5*3  and  a  transposed  form  of  5.13*  the 
following  holds  for  an  arbitrary  positive  numbers  £  o 
and  £  ., 
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*   J)  *         l  6 


In  order  to  obtain  a  non-trivial  bound  we  choose 


5.17  £     =  m    tr* 

3  l     o 

and  hence 

5.18  tfff^t^jJv 


D 


So  that  the  last  equation  may  be  employed  to  bound  5 .10 
we  rewrite  the  previous  equation  as  follows 

5.19     itVs  ^(^jJWA^  JTZ  ^  Jv 

and  then  choose  £ .  so  that  the  coefficients  on  the  right 
are  equal.   We  solve  the  equation 

5.20         2  I  *   +  2MX  £  x  -  m  M2  =  0 

and  obtain 
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5.21 


£   =  k"1   -21^  +^M?  +  8m  M, 


Since  £   is  to  be  a  strictly  positive  number,  the  positive 
square  root  is  chosen  and  thus 


5.22 


£ x  =  2  (-M-l  +  yMJ  +  2m  Mg) 


To  shorten  the  notation  we  let 


5.23 


->       m,  M. 


1*2   L 
and  bound  5»19  with  5>.l8   as   follows 


|  (-M^  +  ^Btj2*  2m1M2)  +  M 


5.24 


Since  a  non-trival  bound  of  the  left  side  of  the  last 
equation  is  desired  we  choose  the  arbitrary  number  c 
in  an  appropriate  manner*. 

I 


$.2$ 


£  =-2-° 
2   M3 


Substitution  of  5.25  into  5»2l}.  yields 

5.26  UUs<£t;ffetJv<$$L'-Js  -  iS 

which  gives  an  estimate  for  the  unknown  boundary  term  of 
$»^>   when  this  term  is  first  bounded  by  5*7. 

To  bound  the  unknown  volume  term  of  5*5  the  following 
cases  are  considered.   For  N  =  2  or  N  =  3  the  results  of 
pages  5o2   and  5&3  °**  [2l]  a**e  applicable  since  from  5»l8 
and  5*26  we  have 


12? 

s.zi  fit  %fh  <  ^  r%JM  i  {$t\Hhf=V 

For  N  >3  a  function     ~h  is   defined  by    (lj..l8)    of    [2lJ    so   that 

5.28  "^F  *   M  x€D 
and 

5.29  BC^)  <  0  x6D, 
Prom  the  divergence  theorem,  it  can  be  shown  that 

5.30  +SS{s&V6Cf)  -  <PX  B&OJjv  < 


<    l  MA  V  +   tLtilH&'SJs} 


k 


"**  -} 

where 
5.31 

M,     =  max  ] 
4          B     L    ^ 

Hence 

5.32 

f  fi>*ff(T)J 

bx^l . 


^IIR^JW^ 


+{j/r-5fr3r*vA  f1^7* 
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The  desired  result  is  now  obtained  by  applying  Schwarz ' s 
inequality  to  the  unknown  volume  Integral  which  we  are 
seeking  to  bound.   We  obtain 

and  note  that  the  first  factor  on  the  right  is  bounded  by„f>.32 
while  the  second  is  computable  due  to  (l|.l)  of  [2l]  . 
Prom  the  above  it  is  clear  that 


5.33  \l^xUrf)AM-Us&)Mvf{SSm)y[KrrfJ<vf 


5-3k  \nt>)-X(p)\, 


*  mk{$  (&h;f%}\  1  qrstQ  h 


where  U  is  given  by  5>»26,  the  second  term  on  the  right  is 
bounded  by  5>«33  a^d 

5.35  i(f>)  =  -3rf>Lels  +  SSrfFjv  +  ifcf>)  ■ 

2.   A  Generalized  Dirichlet  Problem 

In  this  section  a  pointwise  bound  for  the  solution 
of  a  Dirichlet  problem  shall  be  obtained  where  the  given 
differential  operator  is  such  that  the  results  of  f2lj  are 
not  applicable. 

Since  the  problem  and  the  techniques  which  are  to  be 
employed  are  very  similar  to  those  which  were  used  for  the 
Degenerate  Dirichlet  Problem,  the  notation  of  J2l]  shall  be 


129 


altered  so  that  it  more  closely  follows  that  which  was 
previously  developed  in  this  paper* 

We  seek  a  pointwise  bound  for  the  solution  of 


5.36 


J(  w)  =  (a±Jw  .)  ,-qw  =  f  (x)       x£D 
9  -1-  >  J 


fr(  w)  =  (aij* 
w(x)  =  h(x)  x  €  B 


where  D  and  B  are  as  defined  in  the  first  section  of  this 
chapter.   The  solution  function  w  is  assumed  to  be  continu- 
ously differentiable  in  D  and  to  have  a  continuous  second 

derivative  in  the  interior  of  a  finite  number  of  regions  the 

i  1 
sum  of  which  is  D«   The  symmetric  matrix  a  J  is  piecewise 

continuously  differentiable  and  satisfies 

5.37      oso.£$N«'*f.*  SflZf?- 


for  all  real  numbers  (£  -,  ,  .  .  .  ]f„)  where  equality  holds 
on  the  left  if  and  only  if  all  the  numbers  are  zero.   The 
bounded  integrable  function  q  =  q(x)  is  very  similar  to  the 
analogous  function  of  the  Degenerate  Dirichlet  Problem?  that 
is,  q  has  the  form 

V 


5.38         q^^q^lE1 


r=/ 


where 


A 
1=1 


5.39  X  t1  - 1 


and  for  N  2  3 
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Conditions  similar  to  lj..8  a  and  l±.Q   b  may  easily  be  obtained 
for  N  =  1  and  N  =  2.   To  shorten  the  argument,  we  consider 
in  complete  detail  only  the  case  N^3«   As  for  the  degenerate 
problem,  we  define  in  as 


5.1+1         0  $  -I*  £JmT  fa  (N-2)2 

I  L  o      J 


■1 

=  m  <1 


So  that  the  divergence  theorem  and  Schwarz's  inequality  are 
applicable  we  assume  that  f  is  integrable  and  square  integrable 
over  D. 

To  obtain  the  desired  bound  for  the  solution  at  an 
interior  point  p  6  D  we  choose  an  arbitrary  function  $    -  tf (x) 
which  is  such  that  J(  If  )  and  If  (x)  approximate  f  (x)  and  h(x) 
as  appropriate.   We  further  require  that  ^  be  twice  piecewise 
continuously  differentiable  in  D.   We  introduce  the  function 

$.\\2.  f(x)    =  w(x)  -  4>  (x) 

and  compute 

(j(f  )    =  F(x)  xfD 

5.1+3        i    , 

I  Y  (x)  =  H(x)  x6  B. 

From  (6.8)  of  [2l]  we  have 
5.1+1+         W(p)-iP(p) 


^/{^^-^(VJA.if^-ij.H}^ 
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where 

SM  »(f,p)-C^PPfl),J-qr,p 

with  Fp  given  by  (l±»2)    of  the  work  by  Payne  and  Weinberger. 
The  unknown  terms  on  the  right  of  $.l\l\.   are  the  second  and 
fourth.   The  final  term  is  the  first  to  be  estimated.   With 
B  replacing  B(yQ)  and  D  replacing  D(yQ}  we  note  that 
equations  1^.78  to  I4.8I  are  valid  for  the  present  problem. 
Next  we  use 

S.ii6  (Mm  H(*VW«fe  t]£  SScl%  %Jm 

0  B  *  N  qo    0  x    '* 

to  write 


where 


5.L8  m,  =  min  n~     f  n^ 

J  B  K 


5-1*9  R2  ■  min  max    £  (x1-^1) 

x  feD  x  £B    ' 


5.50  M1    =  max    SI  (f1) 


5.51  8^   =  max    {  q2  j 
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and,  as  indicated  above,  /C  is  given  by  Ij. .80.   To  determine 
a  bound  for  the  unknown  term  on  the  right  of  5.47 >  the 
divergence  theorem  is  used  as  follows 

i.&  fjwo h  --  -  y ^tPfyii  +  §*  Uh  -Jfy  <rlJv 

D    J*  V 


where  we  have  use  of  an  appropriate  form  of  lj..lj.O.   Trans- 
position of  5>.f>2  yields 

5.53       f/-*-^)0i\^^ 


8  ft         0 


for  any  positive  numbers  Ci      and  cC  .   Since  we  are  seeking 
a  bound  for  the  integral  of  the  square  of  the  conormal 
derivative  over  B  we  are  led  to  choose  the  arbitrary 
constants  in  the  last  inequality  as  follows 

H2a0 


5.  A         «x  -   (1-m)  j^ 
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and 

-         4%^  ) 


d       k   1  1        a0    IT^a 


Substitution  of  these  results  into  5-^-7  yields 

5.56       f(VO^S  <    13 

s 

where  w  is  an  appropriate  computable  constant  which  is  such 
tnat 

5.^7        13=o 

if 

5. £8       p  =  H  =  0 

For  N  =  2  or  N  =  3  we  use  "the  results  of  pages  5>62 
and  5^3  of  [21]  to  estimate  the  unknown  volume  term  in  %.\\\\<. 
These  results  are  directly  applicable  since  £.53  provides  a 
bound  for 

N 

5.59      JJX^iv, 

If  N>3   then  ^  is   defined  by    (i^.18)    of     [21]    so   that 
5.60  ~T$n       ^       J*(<f )    <  0  x6D 

where 

5.61  j*(*)   =  (.'Jf.)  ,-aq^. 

J  J       9  J- 
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The  divergence  theorem  may  be  used  to  establish  that 

5.62  JT{£F>  Tit)-  \f  T*fF)}  «l  V 


&M&ti{K}k  {iv^Jst 


2 


and,  since  the  right  side  of  this  equation  is  computable, 
it  is  evident  that  the  desired  pointwise  bound  is  available 
using  the  techniques  and  the  results  of  the  paper  by  Payne 
and  Weinberger.  In  the  above  equation  $  is  given  by  5„5>6 
while  'v  is  the  bound  for  the  integral  of  a1 J  f    .   f    . 

jl  9  J 

which  is   obtained   from  5>«53   after  application  of  5>«5>6° 


CHAPTER  VI 
EXAMPLE 

1.   Explicit  Non-Linear  Mixed  Problem 

In  this  section  a  bound  shall  be  obtained  for  the 
solution  of  a  non-linear  boundary  value  problem  which  is  of 
the  type  considered  in  the  first  section  of  Chapter  III. 

We  suppose  that  V  is  a  three  dimensional  domain 
bounded  by  the  hyperplanes  y  =  o  and  y  =  1  and  by  the  surface 
whose  equation  is 


l   ,  l2    22 

6.1         xx  +  x^  - 


The  boundary  value  problem  is  explicitly  stated  as  follows: 

6.2  (x^)6Vr«<W< 


oo 


w(x,o)»o  (*,$eD(p)     j  ^  =  /   (*,<$)  es 


where 


6.3        A>o. 

From  the  above  we  obtain 
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e.k 


*V 


a;    =  o 


pi 

M     -  o 


nt  --o 


T*t   ~  I 


*o     -    I 


4>    =  / 


^ 


By  inspection  one  observes  that  the  differential 
operator  of  6.2  is  slightly  different  from  that  of  3.1. 
The  remarks  which  are  found  at  the  conclusion  of  Chapter 
III  suggest  the  reason  for  this  difference.   The  additional 
term  is  combined  with  the  differential  operator  in  order  to 
indicate  the  improvement  which  follows  from  such  an 
association.   If  this  were  not  done,  then 


6.1j.  a 


"^    1 
1   A2 


and  by  close  inspection  of  the  results  which  are  obtained 
below,  it  is  evident  that  the  sharpness  of  the  error  bound 
is  materially  diminished. 

Following  the  techniques  of  the  earlier  chapters  we 
let 


6.5 


^  =  |r2 


and  then  calculate 
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F  =  (A  +  |  r2)"1 
6.6        ^G  =  «|  r2 

L  =  0 

Even  though  LP    satisfies  the  boundary  conditions  on  S  and 
for  large  A  approximates  the  solution  of  the  differential 
equation,  it  is  evident  that  tf    does  not  closely  coincide 
with  w  throughout  V  since  on  D(o)  the  approximation  is  very 
poor. 

The  set  of  auxiliary  functions  which  were  introduced 
in  the  remarks  accompanying  equation  3»23  may  he  defined  as 
follows: 

i  =  1  and  2 


6.7 

f  1   =  x 

Thus  we 

have 

K  =  1 

6.8 

< 

M1   =  1 

Mo  =  2 

M3 

=  0 

M^ 

=  0 

M^ 

=  0 

In  what  follows  the  solution  of  6.2  is  bounded  at 
the  point  p  =  (0,1).   Prom  3«l6  we  have 

6.9        Iw(o,/)-z(o./)U|j2Uu  JttVJ*(tf„)*lv| 

+  |£jw.  i JJ  Vp  {-f  (xj  ^  Wj  rw) -f^j^Vif)}  J  V  I 

«  -S 
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where 

6.10 

tfp   -    [l+TT  (l-y)J  -1   exp  [-r2  (^(l-y))-1" 

and 

*  * 

6.11 

*  V  -&*-*&-  V  "$l 

Since  If    is  such  that  L  is  zero  on  Ss    we  find  it 
advantageous  to  choose  the  arbitrary  constants  of  3»35>  in 
a  slightly  different  manner  than  in  equations  3«3°'  to  3»42 
Thus  by  setting 


6.12 


oC2   =a> 


and 


6.13 


equation  3»3£  may  be  rewritten  as 

6.14 


fJjt^Av  f  («,,-£  -  $i)  Jtff^v 


V 


V 


9  D[o) 


+  **3  ill(*  +  i^)~^[^(/-^)]i\/ 


Substitution  of   this    into   3»2l|   yields 
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6.15       IJfjs  <  W^fH    f  gfcjlK^**f»W>3 


Due  to  the  form  of  the  above  we  set 

6.16  l|:-'^f^?^ 

and  then  seek  to  determine  the  positive  constants  cC  ^  and 
o( -  so  that  the  right  side  of  6.15  is  minimized.   To  solve 
such  a  minimum  problem  is  extremely  difficult;  hence  in 
practice  we  ordinarily  choose  the  arbitrary  constants  in 
some  admissible  manner  and  then  seek  to  improve  our  result 
by  a  better  choice  of  the  approximation  function  (f   .   For 
example  we  assume  that 


6.17 

A2      16 

and   set 

6.18 

cC1*k 

and 

6.19 

°C,=  <tfj 

oC3=^0-^'^"x)  "  %~  TK   • 


Prom  the  above  and  6.l6  it  follows  that 
6.20         K  =  0. 


ll+O 


Equation  6,15  may  now  be  rewritten  as 

6.21  tifcts  <  7  t      im**o(fr-$,)    =    18 
We  use   6.10   to  obtain 

6.22  (IK^)^:    {W».W 


and  hence 


6-23        ttJV<*s}<t{ir($,)\fc}<1 


«  .3?v  13 


fk 


where  |J  is  given  by  6.21. 

To  bound  the  first  term  on  the  right  of  6.9  we 
introduce  an  auxiliary  function  r   which  Is  such  that  the 
conclusions  of  Theorem  2.2  are  valid.   The  auxiliary 
function  used  in  this  example  Is  not  identical  to  that 
defined  by  2.39.   While  the  original  function  could  be 
used  for  the  present  problem,  advantage  is  made  of  the 
simplicity  of  the  differential  operator  to  obtain  a  function 
which  lends  itself  to  the  integrations  which  are  required 
for  the  determination  of  the  desired  bound.   Let 

6.24        ^fm   (l-y)"£  exp  [-r^U-y))"1]  , 
Clearly 


ll\.l 


6.25  ^0  for  (x,y)  6  V 
and  since 

6.26  T(T)  B  -^  +  — t  +  -^ 


for  (x,y)  €  V,  it  follows  that  the  conclusions  of  Theorem 
2.2  are  valid  as  was  indicated  above.   For  the  present 
differential  operator  equation  3.57  becomes 

6-27        &U.   10  f>T  W($0  -  ^W) 

■ 

s  Do)  r     y  ^ 


Prom  3.62  and  the  above  we  have 

6.28  K  =  K2  ^  (-1  +  ~2)  <0 

where  the  strict  inequality  holds  on  the  right  due  to  6 ,17 
In  what  follows  we  set 

6.29  K2  =  0 
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to  simplify  the  numerical  calculations.   (If  greater 
accuracy  were  desired,  one  would  of  course  take  K2  equal 
to  A~'  -1)-   Prom  equations  3.57  to  3«^5  we  obtain 

6.30  Jt/{^  v  f  -  ^j e=r)]  ^v 

< -6/^1- If]  *-.*/?  =  ^ 


To  use  3*^9  it  is  also  necessary  to  compute 


«  £  tfJo-^/f  *x*p  [-^r/-^))"'J  A\/ 


=  3.aV  fl 


-z 


and 


6.32 


v  r 


=  (v^fU/i)"^p(y(^))"]^ 


=  O.OH/Z 


By  combining  the  foregoing  and  equation  3*69,  it;  follows 
that 
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6-33  \JU,    tff  f  j*(v  )  fa 


<o,iov(/.9ft'l+{f.iii7t-l.l+%i-     /qpr^WTXAJf  J.2V/9  J  /. 


To  bound  the  second  term  on  the  right  of  6.9s  we 
could  use  the  techniques  associated  with  equations  3«&9 
to  3«75?  however,  since 


6.34        J*Ut,)  -  -IT 


P 


we  apply  Schwarz's  Inequality  to  the  expression  under  con- 
sideration and  write 

aj-*/     vfo) 


<  *  { jr/c- wr  j  va^Ym-  iefa  ^/vi 


*i 


+  3,2y  *"*  r.Zn}    /  . 


The  bound  for  the  solution  is  completed  by  calculating 
the  function  X (p)  which  Is  given  by  3»17»   One  easily 
obtains 

6.36  if(o>i)   +-  SS  Mpjs     -  o 

6.37  iT  &  Y*  Is   =  -cog 

D(o)         f 
and 

6.38  -.22/Z/f\<-iI(^F<j/v  *  "  22l2(pf-fe)"f 

Prom  6.9  and  the  above  we  write  as  the  desired  bound  for 
the  solution 

6.39  -Q.38*  {%  +   *W.)(a  A \}^  -a 20V  (/  » j"')^/, gi"' 
f  {+.2/7  f  U*[?  +  ^tO(g-l^3ay^"Z^) 


+  0,0(0   +■  ,ju/z  (a  +  i) 


-a. 


«  w(o0/) 


f-0,ok   +•,££/?.  /T*\ 
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For  example  let  U3  suppose  that 

6.i|0         A  =  1| 

then  one  obtains  from  6.42  the  following  result 

6.^1       ~o.95*w(o,i}  i  +1.09 

When  bounds  for  the  solution  of  more  complicated 
problems  are  desired,  the  techniques  are  exactly  the  same 
as  used  above;  however  the  integrals  which  result  may  be 
quite  difficult  to  evaluate . 
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